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ABSTRACT 


A  finite  difference  time  domain  (FDTD)  technique  for  the  analysis  of  ultrawideband 
(UWB)  antennas  with  rotational  symmetry  has  been  developed.  The  theory  is 
intended  to  treat  various  types  of  biconical  antennas,  and  a  mode  matching  analysis 
of  spherically-capped  bicones  has  been  performed  to  verify  the  correct 
implementation  of  the  FDTD  calculation.  In  addition,  experimental  results  have  been 
obtained  which  are  compared  to  the  theoretical  results.  The  FDTD  software  can  be 
used  for  the  design  of  high-power,  directive  UWB  antennas  with  low  sidelobe  leakage 
for  high-power  UWB  radar  and  electronic  countermeasures  applications. 


RELEASE  LIMITATION 


ocess  additional  to  tb 
USafttUUC  Others 


distribution  is  Ihuitgdto  the  DefeiffeTtonujiun itff  of  Australia, 
’UST berefbvedjp  Cfjief,  ElectftyrdcJdlaj/are  DivisionTESRL. 


Department  of  Defence 
- 4 - 


DEFENCE  science  and  technology  organisation 


DSTO-RR-0019 


Published  by 

DSTO  Electronics  and  Surveillance  Research  Laboratory 
PO  Box  1500 

Salisbury  South  Australia  5108 

Telephone:  (08)  259  5181 
Fax:  (08)  259  5938 

©  Commonwealth  of  Australia  1994 

AR-008-994 

February  1995 


Conditions  of  Release  and  Disposal 

1.  This  document  is  the  property  of  the  Australian  Government;  the  information  it 
contains  is  released  for  defence  purposes  only  and  must  not  be  disseminated 
beyond  the  stated  distribution  without  prior  approval. 

2.  The  document  and  information  it  contains  must  be  handled  in  accordance  with 
security  regulations  applying  in  the  country  of  lodgement,  downgrading 
instructions  must  be  observed  and  delimitation  is  only  with  specific  approval  of 
the  Releasing  Authority  as  given  in  the  Secondary  Distribution  statement. 

3.  This  information  may  be  subject  to  privately  owned  rights. 

4.  The  officer  in  possession  of  this  document  is  responsible  for  its  safe  custody. 
When  no  longer  required  this  document  should  be  destroyed  and  notification  sent 
to:  Senior  Librarian,  Defence  Science  and  Technology  Organisation  Salisbury 
Research  Library. 


DSTO-RR-0019 


A  Finite  Difference  Time  Domain  Calculation  for 
Rotationally  Symmetric  Ultrawideband  Antennas 


EXECUTIVE  SUMMARY 

A  finite  difference  time  domain  (FDTD)  technique  has  been  developed  for  analysing 
rotationally  symmetric  ultrawideband  (UWB)  antennas  of  arbitrary  cross  section. 
This  includes  conventional  biconical  antennas.  The  restriction  to  rotationally 
symmetric  antennas  means  that  a  2  dimensional  calculation  is  possible,  reducing  the 
required  computational  effort  for  realistic  antennas  to  within  practical  limits. 

For  high  power  UWB  applications  such  as  impulse  radar  and  disruption  of  electronic 
systems,  compact,  directive  antennas  with  high  breakdown  voltages  are  required. 
The  FDTD  approach  allows  both  the  antenna  and  its  feed  to  be  accurately  modelled. 
The  technique  can  also  be  used  to  predict  the  characteristics  of  narrowband  antennas 
for  more  conventional  applications. 

The  correct  implementation  of  the  theory  is  verified  by  the  analysis  of  a  9  mm  conical 
D-dot  sensor  used  for  transient  electromagnetics  research  at  DSTO.  Measured  results 
are  presented  which  agree  well  with  the  predictions  of  the  FDTD  analysis.  In 
addition,  a  mode  matching  technique  capable  of  analysing  conventional  biconical 
antennas  is  used  to  check  the  FDTD  calculation.  For  the  conical  D-dot  sensor,  the 
agreement  between  the  FDTD  and  mode  matching  techniques  is  excellent. 


DSTO-RR-0019 


THIS  PAGE  INTENTIONALLY  BLANK 


DSTO-RR-0019 


Authors 


Geoff  Staines 
Electronic  Warfare  Division 

Geoff  Staines  received  the  BSc  (Physics)  degree  from  the 
University  of  Queensland  in  1987,  the  BSc(Hons)  from 
Flinders  University  in  1988,  and  a  PhD  in  plasma  physics 
from  Flinders  University  in  1991.  Since  1992,  he  has  been  a 
research  scientist  in  Advanced  Concepts  Group,  Electronic 
Warfare  Division,  DSTO  Salisbury. 


Sean  Braidwood 

Electronic  Warfare  Division 

Sean  Braidwood  received  the  BSc. (Hons)  from  Adelaide 
University  in  1989  and  a  MSc.  in  atomic  physics  from 
Flinders  University  in  1993.  He  is  currently  working  for 
Advanced  Concepts  Group  in  Electronic  Warfare  Division, 
DSTO  Salisbury. 


DSTO-RR-0019 


THIS  PAGE  INTENHONALLY  BLANK 


DSTO-RR-0019 


Contents 

1.  INTRODUCTION  1 

2.  THE  FINITE  DIFFERENCE  TIME  DOMAIN  (FDTD)  METHOD  2 

2.1  Introduction  2 

2.2  The  Yee  Algorithm  2 

2.2.1  Mesh  Equations  2 

2.2.2  Boundary  Conditions  4 

2. 2. 2.1  Edge  Boundary  Conditions  4 

2.2. 2. 2  Antenna  Feed  7 

2.2. 2. 3  Conducting  Surface  Boundary  Conditions  9 

2.3  The  Contour  Path  (CP)  Approach  10 

2.3.1  Mesh  Equations  10 

2.3.2  Boundary  Conditions  13 

2.3. 2.1  Horizontal  Cell  Deformation  14 

2. 3. 2.1  Vertical  Cell  Deformation  15 

2.5  Far  Field  Calculation  17 

2.6  Energy  Calculation  24 

2.7  Coarse  Outer  Mesh  25 

2.8  Summary  27 

3.  CONICAL  MONOPOLE  D-DOT  SENSOR  ANALYSIS  27 

3.1  Introduction  27 

3.2  Reflected  Voltage  28 

3.3  Transmit  Response  31 

3.4  Receive  Response  34 

3.5  Comparison  of  FDTD  D-dot  Response  with  Experiment  37 

4.  CONCLUSIONS  39 

APPENDIX  A 

SPHERICALLY  CAPPED  CONICAL  MONOPOLE  ANALYSIS  41 

A.l  Spherical  Eigenmodes  42 

A. 2  Aperture  Field  Matching  48 

A. 3  Power  Conservation  53 

A. 4  Small  Conical  Monopole  Approximation  55 

A. 5  Transfer  Functions  57 

A. 6  Time  Domain  Response  60 

A.  7  Summary  62 

A.  8  References  63 


DSTO-RR-0019 


Figure  2.1:  Interleaved  electric  and  magnetic  fields  for  FDTD  grid  cell.  3 

Figure  2.2:  Truncated  mesh  in  cylindrical  geometry  5 

Figure  2.3:  Staircased  conducting  boundaries  9 

Figure  2.4:  Ampere  and  Faraday  contour  geometry  for  Contour  Path  technique  1 1 

Figure  2.5  Florizontal  cell  deformation  near  conducting  surface  1 5 

Figure  2.6:  Vertical  cell  deformation  near  conducting  surface  1 6 

Figure  2.7:  Fine/coarse  mesh  interface  detail  26 

Figure  3.1:  Conical  monopole  geometry  28 

Figure  3.2:  Reflected  voltage  from  ideal,  full  analytic,  and  FDTD 

calculations  for  tp  =  50  ps  30 

Figure  3.3:  Reflected  voltage  from  ideal,  full  analytic,  and  FDTD 

calculations  for  tp  =  100  ps  30 

Figure  3.4:  Reflected  voltage  from  ideal,  full  analytic,  and  FDTD 

calculations  for  tp  =  250  ps  31 

Figure  3.5:  FDTD,  full  and  ideal  analytic,  and  first  order  far  field 

comparison  for  tp  =  50  ps  32 

Figure  3.6:  FDTD,  full  and  ideal  analytic,  and  first  order  far  field 

comparison  for  tp  =  100  ps  33 

Figure  3.7:  FDTD,  full  and  ideal  analytic,  and  first  order  far  field 

comparison  for  tp  =  250  ps  33 

Figure  3.8:  FDTD,  full  and  ideal  analytic,  and  first  order  receive 

voltage  for  tp  =  50  ps  35 

Figure  3.9:  FDTD,  full  and  ideal  analytic,  and  first  order  receive 

voltage  for  tp  =  100  ps  36 

Figure  3.10:  FDTD,  full  and  ideal  analytic,  and  first  order  receive 

voltage  for  tp  =  250  ps  36 

Figure  3.11:  Measured  and  predicted  conical  D-dot  sensor  output  voltage  37 

Figure  3.12:  Reconstructed  and  predicted  incident  field  38 

Figure  A.l:  Conical  monopole  geometry  41 

Figure  A.l:  Transmission  transfer  function  of  matched  conical  monopole  59 

Figure  A.3:  Reception  transfer  function  per  unit  antenna  length 

for  matched  conical  monopole  59 

Figure  A.4:  Radiated  field  from  50  W  conical  monopole  for  various  a/ tw.  61 

Figure  A.5:  Voltage  across  50  W  conical  monopole  receive  antenna 

for  various  a/tw.  62 


DSTO-RR-0019 


1.  Introduction 

High-power  ultxawideband  (UWB)  systems  are  under  consideration  for  applications 
in  radar  and  electronic  countermeasures  [1-4].  Modern  electronic  systems  commonly 
used  for  military  purposes  are  increasingly  vulnerable  to  high-power  emissions  such 
as  UWB  radiation  due  to  the  increasing  use  of  large  scale  integrated  circuits  and  other 
solid  state  components. 

High-power  UWB  radiation  is  generated  using  a  suitable  broadband  antenna  excited 
by  a  voltage  transient.  The  essential  component  of  the  generator  used  to  provide  the 
electrical  transient  is  a  high-power  switch  capable  of  producing  subnanosecond 
output  pulse  risetimes.  Photoconductive  solid  state  switches  [1,  3,  5,  6]  or  spark  gap 
devices  [1,  7,  8]  are  suitable  for  this  purpose.  While  high-power  UWB  generators  are 
both  available  and  affordable,  one  of  the  most  critical  hardware  items  for  a  practical 
UWB  system  is  the  antenna.  There  are  several  important  requirements  for  a  useable 
antenna: 


(i)  High  breakdown  voltage  (>  100  kV)  in  order  to  radiate  sufficiently 
large  voltage  transients  to  detect  (for  radar)  or  disrupt  (for  countermeasures)  targets 
at  useful  ranges.  This  apparently  ambitious  requirement  is  mitigated  somewhat 
since  the  voltage  peaks  across  the  antenna  for  approximately  a  nanosecond  or  less, 
and  so  the  pulse  breakdown  voltage  is  up  to  an  order  of  magnitude  higher  than  the 
DC  breakdown  voltage. 

(ii)  High  directivity  in  order  to  maximise  the  power  incident  on  the  target, 
while  avoiding  irradiation  of  nearby  friendly  or  neutral  targets.  The  directivity  is 
ultimately  limited  by  the  maximum  permissible  size  of  the  antenna  aperture  as 
dictated  by  host  platform  constraints. 

(iii)  Low  sidelobes  in  order  to  minimise  the  possibility  of  damage  to  or 
disruption  of  the  host  platform.  While  hardening  of  the  host  platform  is  possible,  the 
cost  and  effort  required  is  minimised  by  the  use  of  a  low-sidelobe  antenna.  This  goal 
is  the  most  difficult  to  achieve  since  the  antenna  must  cover  a  very  large  bandwidth 
(typically  greater  than  25%  of  the  centre  frequency),  while  being  non-dispersive  (no 
distortion  of  pulse  shape). 

Since  fully  three-dimensional  electromagnetic  field  calculations  require  a  substantial 
computational  effort,  effectively  two-dimensional  rotationally  symmetric  antennas 
such  as  the  bicone  were  chosen  for  detailed  study.  All  of  the  previously  mentioned 
aspects  of  high  power  UWB  antenna  design  can  be  investigated  using  a 
straightforward  finite  difference  time  domain  (FDTD)  algorithm  requiring  a  modest 
computational  effort.  This  report  describes  the  development  of  this  FDTD  technique, 
together  with  initial  results  relating  to  the  analysis  of  small  spherically-capped  conical 
monopole  D-dot  sensors.  These  results  include  a  comparison  of  results  obtained 
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from  the  FDTD  analysis  and  an  analytic  mode  matching  technique  for  these  sensors. 
Experimental  results  are  included  to  show  that  the  FDTD  technique  successfully 
predicts  the  correct  dependence  of  the  sensor  output  voltage  on  the  angle  of  incidence 
of  the  applied  field.  A  detailed  derivation  of  the  analytic  technique  is  presented  in 
Appendix  A.l. 


2.  The  Finite  Difference  Time  Domain  (FDTD) 

Method 


2.1  Introduction 

There  are  currently  two  categories  of  time  domain  numerical  techniques  commonly 
used  for  analysing  UWB  antennas.  The  first  type  are  the  integral  equation  methods 
such  as  the  electric  field  integral  equation  (EFIE)  method  developed  by  Rao  [9,  10]. 
Conducting  surfaces  are  modelled  by  planar  triangular  or  rectangular  patches,  and 
the  method  of  moments  is  applied  in  the  time  domain  to  calculate  the  current 
densities  over  these  surfaces.  The  second  type  consists  of  the  finite  difference  time 
domain  (FDTD)  method  first  proposed  by  Yee  [11].  This  technique  is  used  to 
determine  the  fields  in  a  volume  surrounded  by  a  surface  on  which  a  boundary 
condition  is  imposed.  Maxwell's  equations  are  discretised  spatially  and  temporally, 
and  the  fields  in  the  volume  are  solved  on  a  time-marching  finite  difference  mesh. 
The  FDTD  method  has  the  advantages  of  simplicity  and  generality  over  integral 
equation  methods,  although  at  the  expense  of  a  larger  computational  effort,  since  the 
electromagnetic  field  in  a  three-dimensional  volume  is  being  calculated  as  opposed  to 
the  current  density  over  two-dimensional  surfaces.  The  memory  requirement  is 
alleviated  by  considering  only  two  dimensional  azimuthally  symmetric  antenna 
structures. 

2.2  The  Yee  Algorithm 

2.2.1  Mesh  Equations 

Since  the  antennas  are  rotationally  symmetric,  the  analysis  can  be  reduced  to  a  two- 
dimensional  problem  by  adopting  cylindrical  coordinates  [12-14],  Since  the  excitation 
for  the  rotationally-symmetric  antenna  is  a  TEM  mode,  only  TM  modes  will  be 
excited.  In  cylindrical  geometry,  these  modes  have  field  components  'Er,  Ez,  and 
only.  Maxwell's  time  dependent  curl  equations  relating  these  quantities  are 


dEr  dEz 

9H* 

dz  "  dr 

■  dt 

(1) 

dH<j, 

3z 

dEr 

~  E  dt 

(2) 
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1 3frH+)  _  3Ej 

r  dr  £  3f  W 

where  e  is  the  permittivity  at  the  field  point,  and  po  is  the  permeability. 

The  field  quantities  are  discretised  in  the  spatial  and  time  coordinates  for 
representation  on  a  finite  difference  mesh.  For  spatial  increments  Ar,  and  Az,  and  time 
increment  Af,  the  notation  for  field  quantity  Q  at  the  point  (rf,  Zj)  at  time  tn  is 

Q  {hj)  =  Q(ri,  zp  t „) 

=  Q(iAr,  ;AZ,  nAf) 


e”  (/,/'+ 1/2) 


Figure  2.1:  Interleaved  electric  and  magnetic  fields  for  FDTD  grid  cell. 


To  implement  the  Yee  algorithm,  the  electric  and  magnetic  field  components  are 
evaluated  at  interleaved  spatial  grid  points  and  time  steps  [11],  as  shown  in 
Figure  2.1.  The  derivatives  in  Equations  1-3  are  replaced  by  the  second-order 
accurate  centred  difference  approximations,  so  that  Maxwell's  equations  become 


n+l/2 

H 


{hi) 


+  PJ5  [e"(M/2,/)  -  <(i-l/2,/)]  -  jjjs;  [e"(;,;+1/2) 

-  e”(i',/-1/2)]  (4) 


n+1  n  Af  n+l/2  n+l/2 

Er  (i,j+ 1/2)  =  Er(i,;+l/2)  -  ^  [H0  (/,;+!)  -  (i,j)] 


(5) 
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n+ 1  n  Af  n+ 1/2  n+1/2 

Ez  (i+l/2,j)  =  Er(z+l/2,/)  + - [r!+1H^  (i+l,j)  -  [rfH0  (*,/)]  (6) 

ri+  l/2eAr 


The  field  points  HQ  (z, /),  Er  (i,j+ 1/2),  and  Ez  (z+1/2,/)  are  assigned  to  the  grid 
cell  (z,/)  centred  at  (ri,zj) 

For  convergence,  the  mesh  increments  Ar,  A2,  and  Af  must  satisfy  the 
Courant-Friedrichs-Lewy  condition 


c  Af 


2  2,2  2 
Ar  Az/(Ar  +  Az) 


where  c  =  1/ Vepo  is  the  velocity  of  light. 

2.2.2  Boundary  Conditions 
2. 2. 2.1  Edge  Boundary  Conditions 

Practical  considerations  require  that  the  mesh  must  be  finite  in  its  spatial  dimensions 
as  shown  in  Figure  2.2.  From  Equations  4-6,  it  is  apparent  that  to  simulate  an 
unbounded  space,  a  boundary  condition  is  required  at  the  edges  of  the  mesh  to 
specify  the  tangential  component  of  the  electric  field.  At  r  =  0  (i  =  0),  symmetry 
about  the  axis  requires  that 


Therefore 

E"(l/2,;)  =  Ez(-l/2,;') 

and  so 

n+1/2  n- 1/2  Af  n  n 

(0,/)  =  (0/  j)  -  [Er(0,7+l/2)  -  Er(0,;-l/2)]  (7) 

Hence  Ez  (1/2,7)  *s  not  required  explicitly.  Including  the  tangential  electric  field  at  the 

other  three  sides  of  the  mesh  is  not  so  straightforward,  and  must  be  determined  using 
absorbing  boundary  conditions  in  order  to  simulate  an  unbounded  medium. 
Merewether  [14]  used  the  radiation  condition  for  the  electric  field  at  a  sufficiently 
large  distance,  R,  from  the  centre  of  the  antenna. 
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£ 


Figure  2.2:  Truncated  mesh  in  cylindrical  geometry 


_  f {t  -  R/c) 

E  "  R 

where  f  is  some  causal  function,  ie.,  f(x)  =  0  for  x  <  0.  This  corresponds  to  an 
assumption  that  at  a  reasonable  distance  from  the  antenna  the  fields  can  be 
represented  in  terms  of  an  outward  travelling  spherical  wavefront.  Consider  Er  just 
beyond  the  upper  axial  limit  of  the  mesh  (j  =  /,  z  =  zmax) 

R(i, /+1/2)  Ef1(i,/+l/2)  =  fr(tn+l  -R(i,J+l/2)/c) 

=  fr(f”  -  R(i,/-l/2)/ c  +  Af8) 


where 
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c  At 


;  0  <  1 


and 


mi) 


■V 


2  2 

n  +  (zj  -  zo) 


with  Zg  the  axial  position  of  the  geometric  centre  of  the  antenna. 

Similarly,  the  field  just  inside  the  upper  mesh  boundary  (j  =  J ,  z  =  zmax)  is  given  by 

R(i,/-l/2)  E"0,J-l/2)  =  fr(f”-R(z,/-l/2)/c) 

tt+ 1 

Therefore,  it  is  possible  to  obtain  Er  (/,/+ 1/2)  using  parabolic  interpolation  between 
the  elements  of  the  time  sequence  of  e”(z,/-1/2)  as  follows 

n+ 1.  .  .  K(z,/+l/2) .  n- 1, 


Er  *(i,/+l/2)  =  2R^7-l/2)  W0"1)  Er  *(J''/-V2) 


n+ 1. 


+  2(1-02)  Er(z',/-l/2)  +  0(0+1)  Er  (zJ-1/2)] 


(8) 


n+ 1. 


(9) 


Similarly,  Er  (z',-1/2)  at  the  lower  axial  limit  of  the  mesh  (/  =  0)  is 

E;!+1(z,-1/2)  =  2R&1/2)  I0^0-1)  E r*1(U/2)  +  2(1-02)  1/2) 

+  0(0+1)  e”+1(z',1/2)] 

n+1 

Ez  (1+1/ 2,/)  at  the  radial  limit  of  the  mesh  (z  =  I)  is  calculated  in  a  similar  fashion 
starting  from 

R(f+l/2,/)  e”+1(/+1/2,7)  =  f2(f”+1-R(f+l/2,;)/c) 

=  fz(fn-R(/+l/2,;)/c  +  Af0) 


where 


R(/+l/2,;)  -  R(/-l/2,/) 


;0<1 
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Similarly 


R(M/2,;)  e"(M/2,;)  =  fz(f"  -  R(M/2,/)/C) 

/z+1 

Therefore,  it  is  possible  to  calculate  Ez  (1+1/ 2, j)  at  the  side  of  the  mesh  using 

E”+Vl/2,/)  =  2R(!ll2/2,i)  [e(6'1}  E"’Vl/2d)  +  2(1'e2)  eJ(M/2,/) 

+  0(9+1)  EfV-1/2,/)]  (10) 

Equations  7-10  are  sufficient  to  determine  the  required  tangential  electric  fields  on  the 
mesh  boundary,  provided  that  R  is  sufficiently  large  that  the  fields  at  the  mesh 
boundary  can  be  represented  by  spherical  waves  propagating  outward  from  the 
centre  of  the  antenna.  Since  R  must  be  finite  with  a  finite  mesh  size.  Equations  7-10 
constitute  approximately  absorbing  boundary  conditions. 

2.2. 2.2  Antenna  Feed 

Power  is  fed  into  the  antenna  through  a  coaxial  line  extending  directly  downward 
from  the  antenna  feed  point.  The  line  dimensions  are  such  that  it  can  only  support  a 
TEM  mode,  which  means  that  an  exact  absorbing  boundary  condition  can  be 
constructed  inside  the  line.  The  impedance  of  the  line  is 

Zline  = 

and  is  usually  chosen  to  be  50  Q.  The  TEM  mode  electric  field  in  the  line  can  be 
written  as 


Er  =  h(f  +  z/c)  +  g(f  -  z/c )  E,),  =  Ez  =  0 

where  g(f  -  z/  c)  represents  the  TEM  wave  travelling  toward  the  antenna  feed,  and 
h (t  +  z/  c)  represents  the  TEM  wave  reflected  from  the  feed  due  to  imperfect  matching 
of  the  line  to  the  antenna.  This  mismatch  will  also  generate  higher-order  evanescent 
modes,  so  the  coaxial  line  is  extended  a  sufficient  distance  from  the  antenna  feed 
point  such  that  these  evanescent  modes  have  decayed  to  a  negligible  level  before  the 
boundary  condition  is  imposed.  A  distance  of  3 (b-a)  is  ample  for  this  purpose  [12]. 
If  the  feed  line  is  terminated  with  a  match  at  j  =  /term  corresponding  to  z  =  Zferm,  the 
electric  field  inside  the  line  is  given  by 

n+l 

Er  (*// term-1  /  2)  =  h(fn+1)  +  g(f„+1) 


7 


DSTO-RR-0019 


=  h(fB  +  A z/c-  Af0)  +  g (tn  -  A z/c-  Af0) 


So  that 


n+ 1 

Ep  (*//term"l/2)  -  g(f„  -  A2/c  -  A^0)  -  h (tn  +  A2/c  -  Af0) 


n+ 1 

—  ^  (I'//term_l/2) 


where 


0 


and  %(z,f)  represents  the  reflected  TEM  wave.  Proceeding  in  the  same  fashion  as  for 
the  absorbing  boundary  conditions  at  the  extremities  of  the  mesh,  the  reflected  field  at 
the  termination  is 


qn+Vd'term-l/2)  =  ^0(0-1)  ^Vdterm+1/2)  +  2(1-02)  ^(i,yterm+ 1/2) 


+  0(0+1)^"  1(i,/term+l/2)] 


(11) 


To  calculate  ^  (i,j),  the  incident  field  in  the  feed  line  is  subtracted  from  the  total 

field.  The  incident  field  in  the  coaxial  line  is  given  by 


me. 


E r  (0 


Vinc(Q 
In  (b/a)r 


where  Vinc(0  ^  the  user  specified  voltage  applied  to  the  antenna  feed.  A  Gaussian 
pulse  is  commonly  used,  with 


Vinc(0 


r-f2i 


=  V0  exd 


where  xp  =  0.4246  tw,  where  tw  is  the  FWHM  pulse  width. 

Since  only  the  TEM  mode  is  supported  in  the  feed,  (11)  constitutes  an  exactly 
absorbing  boundary  condition.  The  reflected  voltage  on  the  feed  transmission  line  is 
calculated  by  averaging 


Vrefl(0  =  ^(z,;'term-l/2)  ln(b/a)  r 
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over  the  radius  of  the  line.  The  reflected  voltage  constitutes  part  of  the  output  from 
the  program.  It  is  especially  useful  for  identifying  the  nature  of  reflections  from 
various  parts  of  the  antenna  and  for  calculating  antenna  efficiency  for  a  given 
excitation.  The  efficiency  for  a  time  domain  antenna  is  the  ratio  of  the  transmitted  to 
the  excitation  energy  [15],  and  is  calculated  from 

oo  oo 

jVinc(i) 2  d£- jvrefi(02  dt 


/Vine  (t)2dt 
0 

2.2. 2. 3  Conducting  Surface  Boundary  Conditions 

The  implementation  of  conducting  surface  boundary  conditions  using  the  Yee 
formalism  is  particularly  straightforward.  The  surface  is  approximated  by  a 
"staircase"  as  shown  in  Figure  2.3.  Grid  cells  are  filled  if  their  centres  lie  within  the 
conductor.  Er  is  set  to  zero  at  horizontal  interfaces  between  filled  and  unfilled  cells, 
and  Ez  is  set  to  zero  at  vertical  interfaces. 


Ez  = 


<>  +  <> 

I  I 

.  J - • - 1  - 

I  I 

I  I 

o  +  o 

i  i 

- ,  •"  - 

i 

o 


-•  - 
+ 

+ 


Conductor  boundary 


Conductor  boundary 


Figure  2.3:  Staircased  conducting  boundaries 
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Although  simple,  this  technique  is  known  to  introduce  errors  arising  from  the 
staircasing  approximation  when  applied  to  surfaces  which  do  not  coincide  with  grid 
cell  boundaries  [16].  These  errors  are  significant  near  regions  of  high  field 
concentration  such  as  antenna  feed  points.  To  more  accurately  model  the  arbitrarily 
shaped  conducting  boundaries  of  the  antenna,  the  contour  path  formalism  [17] 
described  in  Section  2.3  was  adopted. 

2.3  The  Contour  Path  (CP)  Approach 

2.3.1  Mesh  Equations 

The  contour  path  (CP)  algorithm  described  by  Jurgens  et  al.,  [17]  is  based  on 
Faraday's  and  Ampere's  laws  rather  than  Maxwell's  equations 


fc  Ed'  = 

(12) 

(13) 

where  the  contour  C  encloses  the  surface  S.  The  contours  used  for  Equations  12  and 
13  intersect  each  other's  enclosed  surface  in  a  manner  resembling  the  links  of  a  chain, 
as  shown  in  Figure  2.4. 

Equations  12  and  13  are  implemented  by  defining  the  electric  and  magnetic  field 
values  at  interleaved  grids  in  the  r-z  plane  as  discussed  in  Section  2.2.  As  for  the  Yee 
method,  a  centred  finite  difference  approximation  is  used  for  the  time  derivatives.  To 

evaluate  the  surface  integral  in  Equation  12,  assume  that  is  constant  over  a  grid 
patch  and  equal  to  the  value  at  the  centre  of  the  patch  as  shown  in  Figure  2.4.  In 
addition,  Er  and  Ez  are  assumed  to  have  no  variation  along  the  individual  sides  of  the 
patch.  For  the  cell  centred  at  (ru  zj) 


1 

‘at 


Zi+V 2 
r  r  1+1/2 


a 


r 

./ 


H  dr  dz 


J 

zi-i/2 


n-1/2 


Po  n+ 1/2 

A7(h0  (w) 


n- 1/2 

H*  (i,j))Acen 
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where  Aceu  is  the  area  of  the  grid  cell  in  the  r-z  plane.  From  Figure  2.4,  the  contour 
integral  is  calculated  using 

|  E M  =  E"(*+l/2 ,j)  Ifight  -  e”(z-1/ 2,j)  heft 

+  Er  (1,7-1/ 2)  l [ow  -  Er(i,7+l/2)  lUp 

where  lright  and  l[e^  are  lengths  of  the  left  and  right  sides  of  the  cell,  and  l[aw  and  lUp 
are  the  lengths  of  the  upper  and  lower  sides  of  the  cell.  Therefore 

ti+l/2  n-\/2  At  n  n 

H,  ( i.j r)  -  H,  (i,J)  + (Ez(i+l/2,;)  1,,^,  -  EI(,-l/2,/)  lufi 

+  E"(i,/-l/2)  l!om  -  E°(i,;+l/2)  (14) 

Q 


Figure  2.4:  Ampere  and  Faraday  contour  geometry  for  Contour  Path  technique 


For  a  simple  rectangular  grid  cell,  Aceu  -  ArAz,  lup  =  kow  =  Ar,  and  keft  =  lright  =  A2. 
Equation  14  then  reduces  to  Equation  4  as  derived  from  the  Yee  approach. 

Since  there  is  azimuthal  symmetry,  E”(i,;+l/2)  is  found  from  Equation  13  by 
integrating  over  and  around  the  Ampere  loop  shown  in  Figure  2.4.  Note  that  only  FI,), 
contributes  to  the  line  integral  since  the  contributions  from  the  segment  connecting  the 
two  loops  cancel  each  other.  It  follows  that 
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J1ee- 


is  -  at 


JeEr(z,;+l/2)  (rtd<j))dz 


=  g  (e”+1(z',/+1/2)  -  e”(/,7+1/2)  2m 

diel  diel 

(er(z,;+l)  \lv  +  er(i,j)  \ow) 

diel  diel 

where  and  are  the  distances  from  a  possible  dielectric  interface  to  points 

(rj,  Zj+i)  and  (ruZj),  respectively,  and  zr(i,  j)  is  the  relative  dielectric  constant  at  mesh 
point  ( ri,Zj ).  Note  that  only  a  horizontal  dielectric  interface  is  allowed  at  points 

(n,  Zj+1/2)  where  Er  is  calculated.  The  contour  integral  from  Equation  8  is 

f  n+ 1/2  n+l/2 

jc  H.dZ  =  -2nn(H0  -  H0  (i,j)) 

Therefore 

n+l/2  n+l/2 

n+ 1  „  AfH0  (W+l)-H#  (i,/) 

Er  (i',y+l/2)  =  Er(i,;+l/2)  -  ^ (15) 

£r(N+l)  kip  +  MW)  kow 

diel  diel 

For  the  simple  case  of  uniform  dielectric  constant  and  li0lv+  lUp  =  Az,  Equation  15 
reduces  to  Equation  5  as  for  the  Yee  formalism. 

E  (z,/+l/2)  was  calculated  in  a  similar  fashion.  The  surface  integral  is 


llseE‘ 


dt  JeEz0'+l/2d)  (rid<p)dr 


En  n+ 1  n 

g(Ez  (i+l/2,j)-Ez(i+l/2,j)2K 
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diel  diel 

(Er(i+Xj)  ri+1lright+  n  er(i,j)  lleft) 


where  and  are  the  lengths  from  the  dielectric  interface  to  points  (rI+i,  zj)  and 
(ri,Zj),  respectively.  Note  that  only  a  vertical  dielectric  interface  is  allowed  at  points 

(ri+\/2,  zj)  where  Ez  is  calculated.  The  corresponding  contour  integral  is 


n+l/2  n+l/2 

=  2n(ri+1  (/+!;)  -  rf  (i,j)) 


(  Therefore 

f 


Enz+\i+l/2,j)  =  Enz(i+l/2j) 


n+l/2  n+l/2 

At  ri+l  H0  (i+V)  -  ri  H<p  (W) 

£0  diel  diel 

Er(i+l,j)  ri+11right+  Er  n  (i,j)  lleft 


(16) 


For  the  simple  case  of  uniform  dielectric  constant  and  =  Ar,  Equation  16 

reduces  to  Equation  6  from  the  Yee  formalism. 

It  is  clear  that  for  simple  rectangular  grid  cells,  the  mesh  equations  derived  using  the 
contour  path  approach  are  exactly  the  same  as  those  derived  using  the  Yee  method. 
However,  it  is  near  conducting  boundaries  that  the  contour  path  method  is  most 
useful.  If  one  of  the  sides  of  the  cell  coincides  with  the  conducting  surface,  then  the 
tangential  component  of  the  electric  field  along  that  side  disappears.  It  does  not 
matter  that  the  cell  is  no  longer  rectangular  as  the  shape  of  the  cell  is  not 
fundamentally  constrained  in  the  contour  path  approach  as  it  is  in  the  Yee  formalism. 
A  thorough  discussion  of  cell  deformation  near  conducting  surfaces  is  presented  in 
the  next  section. 

2.3.2  Boundary  Conditions 

The  boundary  conditions  at  the  extremities  of  the  mesh  which  were  derived  for  the 
Yee  formalism  are  equally  valid  for  the  contour  path  method,  so  this  section  will  deal 
exclusively  with  boundary  conditions  near  conductors.  For  the  grid  cells  adjacent  to 
conducting  surfaces,  the  Faraday  contours  used  for  the  integration  in  Equation  12  are 
deformed  so  as  to  conform  to  the  surface  of  the  conductor.  Ampere  contours  are  not 

deformed,  and  Er  and  Ez  values  which  require  using  Ampere  contours  which  cross 

media  boundaries  are  not  used.  For  this  reason,  only  the  calculation  of  values  is 

affected  by  the  presence  of  a  conductor.  The  calculation  of  Er  and  Ez  values  is 
essentially  identical  to  the  Yee  method. 
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2.3. 2.1  Horizontal  Cell  Deformation 

For  the  case  where  a  conducting  boundary  lies  to  the  left  of  a  grid  cell,  and  the  slope 
of  the  conductor  is  greater  than  Az/Ar  with  respect  to  the  horizontal  as  shown  in 

Figure  2.5,  then  is  calculated  from 


H 


n+l/2 

<P 


(w) 


n-1/2  Af 

H*  <''>>*  Mo 

(Ez(i+l/2,;)  l right  +  Er  (i,j-l/2)  Z/ow  -  Er(z,j+l/2)  lup) 


If  the  centre  of  the  cell  immediately  above  located  at  (r;-,  Zj+\)  is  inside  the  conductor, 

then  the  Ampere  contour  used  to  calculate  Er(z,j+l/2)  is  discarded  and  the  nearest 
neighbour  value  on  the  right  is  used,  so  that 

n+l/2  n-1/2 

H0  (i,j)  =  H0  (z ,/) 

+  Po  Aceu  (Ez(Z+l/2,j)  l right  +  E r (z,/-l / 2)  liow  -  Er(z+l,/+l/2)  lUp ) 

Similarly,  if  the  centre  of  the  cell  immediately  below  located  at  (r,-,  Zj. i)  is  inside  the 
conductor,  then 


n+l/2  n-1/2 

(i,j)  =  H,  (Z,/) 

Af 


+  po  Acell  (Ez(z+1/2d)  bright  +  Er  (z+l,/-l/2)  kow  -  Er(z,/+l/2)  lup) 


If  the  conducting  surface  is  to  the  right  of  the  cell,  the  treatment  is  similar  except  that 
the  electric  field  on  the  right  side  of  the  cell  is  set  to  zero.  Nearest  neighbour  values  on 

the  left  are  used  to  interpolate  Er(z,/+l/2)  and  Er(z,;-l/2)  where  Ampere  contours  are 
not  available.  For  this  instance 

„»+l/2 . n-\/2 

C l’l ) 

(e"(z-1/2,;)  lkft  +  e"(z,;-1/2)  llow  -  Er(i,j +1/2)  lup) 


^  (W)  =  ^ 

At 


V-oAcell 
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°  Er 

•  e2 


Figure  2.5  Horizontal  cell  deformation  near  conducting  surface 


2. 3. 2.1  Vertical  Cell  Deformation 

If  the  slope  of  the  conducting  boundary  is  less  than  Az/  Ar  relative  to  the  horizontal,  it 
would  become  possible  for  the  vertical  cell  sides  to  be  truncated,  complicating  the 
situation.  In  addition,  inaccuracy  could  be  introduced  if  the  electric  field  values  are 
held  constant  over  unnecessarily  long  upper  and  lower  cell  sides.  For  these  reasons,  it 
is  better  to  deform  the  grid  cells  vertically  in  this  instance.  For  the  case  where  a 
conducting  boundary  lies  just  above  a  grid  cell  as  shown  in  Figure  2.6,  then  is 
calculated  from 


n+1/2  n-1/2 

H*  (/,;)  =  h,  (i,j) 


At 


+  po  A ceii  (Ez0+l/2,/)  Iright  ~  Ez(z-l/2,/)  lleft  +  Er(i,j-l/2)  Z/ow) 

If  the  centre  of  the  cell  immediately  to  the  right  at  (r,+ zf)  is  inside  the  conductor, 

then  the  Ampere  contour  used  to  calculate  Ez(z+l/2,;)  is  discarded  and  the  nearest 
neighbour  value  from  below  is  used,  so  that 

n+1/2  n-1/2 

H*  (w)  =  H.  (i,j) 


A  t 


j^O  A  cell  bright  ~  E  z(f-l  /  2,j)  Ity \  +  Er(z,/-l/2)  llow) 
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Similarly,  if  the  centre  of  the  cell  immediately  to  the  left  at  (r;_i,  zf)  is  inside  the 
conductor,  then 


n+l/2 

H,  (i,j) 

At 


n-1/2 

h*  di) 


+ 


n . 


hO  A  ceu 


(E2(i+l/2,j)  lrigjit  -  Ez(z-1/ 2,y-l)  +  Er(i+l,y-l/2)  limv ) 


If  the  conducting  surface  is  below  the  grid  cell  the  treatment  is  similar,  except  that  the 
electric  field  on  the  bottom  of  the  cell  is  set  to  zero.  Nearest  neighbour  values  from 

above  are  used  to  interpolate  Ez(z+l/2j)  and  Ez(z-l/2,/)  where  Ampere  contours  are 
not  available.  For  this  situation 


n+l/2 

h,  (hi) 


At 

hO  A  cell 


(hi) 


(Ez(*+1/2'7)  bright  -  Ez(z+l/2,j)  //eyf  -  Er(z,;+l/2)  lup) 


Conductor  boundary 


°  Er 

•  Ez 


Figure  2.6:  Vertical  cell  deformation  near  conducting  surface 
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2.5  Far  Field  Calculation 

The  technique  developed  in  [18]  was  used  to  calculate  the  electric  field  in  the  far  zone 
from  the  fields  close  to  the  object  which  are  determined  using  the  FDTD  method. 
From  [19],  the  radiation  vectors  N  and  L  over  a  closed  surface  S'  with  local  surface 
normal  n  are  defined  in  the  frequency  domain  to  be 


N(k)  = 


L (k)  =  -  [  n  x  E(R',  k)  exp(-;7cR'.R)  dS' 
s"' 

where  k  =  2 nf/c,  R  is  the  vector  from  the  centre  of  the  antenna  to  the  far-field  point 
and  R1  is  the  vector  to  the  point  on  the  surface  S'  enclosing  the  antenna.  Using 
spherical  coordinates  with  the  origin  located  at  the  centre  of  the  antenna,  the  only 
electric  field  components  which  do  not  decrease  faster  than  1/R  are 

e8((t)  -  Coiyt)  -  fVq) 

E,<*)  -  -  Co  He(fc)  -  j  (-  CoN,  +  Lg) 

where  C,o  is  the  impedance  of  free  space.  Following  [18],  define 

Taking  the  inverse  Fourier  transform  of  W  and  U  leads  to 

■  JrcW  jAxH(RI'  n  dS' 


■w  -  -ikie  jAxE^Vds' 


|  n  x  H(R',  k)  exp(-/fcR'.R)  dS' 
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where  V  =  t  -  [  R  -  R'  |  /c  is  the  retarded  time.  The  time-domain  electric  field  in  the  far 
zone  is  calculated  from 

EeOM  =  -(Cowg  +  u^ 

=  ik  ([-  ^  [A  x  H(R'  f)]«  +  if  [fi  x  E(R’,  f')]J  dS‘  (17) 

S' 

E»(R-  0  =  -  S0"'£,  *  u0 

=  5S  j[-  it  ["  x  H(R'  %  -  Je  [n  *  E(RI-  <■)]„]  iS'  (18) 

S' 

Equations  17  and  18  can  be  combined  and  expressed  as 

E<R'  f>  -  is 


J^T[nx  H(R’,  f’)]  dS'+^  j R  x  ~  [A  x  E(R',  f)]  dS' 
S'  S' 


ll(|  f  A  A  3 

=  4^  J  RxRx^[llxH(R',f')]rfS' 

S' 

-h  Jrx^[AxE(R VO]rfS']  (19) 

S' 

Equation  19  is  identical  to  Equation  7  presented  in  [12].  For  convenience,  a 
cylindrical  surface  was  used  to  enclose  the  antenna.  The  side  of  the  cylinder  extended 
from  the  ground  plane  at  j  =  jgp  up  to  j  =  jtop.  The  radius  of  the  cylinder  was  rS/de- 
There  are  contributions  to  the  surface  integral  in  Equation  17  and  18  from  both  the 
side  and  top  of  the  cylindrical  surface.  Converting  from  cylindrical  to  cartesian 
coordinates,  for  the  side  of  the  surface 
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nxH(R ',?)  =  H  (R',  f)  z  and  nxE(R't')  =  -  EZ(R',  0  0 

=  E2(R',  t ')  [sin  0'  x  -  cos  0'  y] 
and  from  the  top  of  the  surface 

n  x  H(R’,  f)  =  -H^(R'  f)r  and  nxE(R'  f)  =  EZ(R',  f1)  0 


=  H^(R',  V)  [  -  cos  0'  x  -  sin  0'  y] 


=  Er(R',  £')  [  -  sin  0'  x  +  cos  0'  y] 


At  the  far  zone  point  where  0  is  set  arbitrarily  to  zero  since  ^  =  0,  then  x,  y,  and  z 

A  A  A 

can  be  expressed  in  terms  of  R,  0,  and  0  at  the  far  zone  point  (R,  0,  0)  as  follows 

x  =  R  sin  6  +  0  cos  6 

a  1 

y  =  0 

z  =  R  cos  0-0  sin  6 
So  at  the  side  of  the  cylindrical  surface 

n  x  H(R';  t')  =  H0(R',  f)  [cos  9  R  -  sin  9  0] 

and  hence 


RxRxnx  H(R',  t')  =  H0(R',  f)  sin  9  0 

Also 

n  x  E(R',  f)  =  E2(R',  £')  [sin  9  sin  0'  R  +  cos  9  sin  0'  0  -  cos  0'  0] 
which  leads  to 

Rxnx  E(R',  £')  =  EZ(R',  £')  [cos  0'  0  +  cos  9  sin  0'  0] 
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The  contribution  to  the  far  zone  electric  field  from  the  side  of  the  cylindrical  surface  is 


side  (In  C  a  a  3  . 

E  (R'0  =  4kR[  J  RxRx^[HxH(R'f)]rfS' 


side 


'Vo  |Rx^[nxE(R',f')]dS'] 

side 


2  K 

r 


ide 

4  kR 


ztop 


sin  6  J  ^rside,  z' ,  f)  dz 1  d<j>'  6 


cgP 


1 


2  n 

P  ztop 


j  ^  E2(R',  f)  dz'  [cos  <f>'  0  +  cos  9  sin  tp'  <j>]  d(f>'  ] 


cgP 


M-0 rside  r  . 

"4^[sm0 


2k 

(*  ztop 

d_ 

dt ' 


H  irside,  z\  f)  dz'  dip'  0 


cgP 


1 

To 


2k 

ztop 

r  a 

I  ^E2(R',  t ')  dz'  cos  d0']  0 
Zo 


(20) 


since  E  (R1,  £')  is  an  even  function  of  <p',  as  can  be  seen  from 


f  =  f- 


R-R' 
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=  t 


sin  9  -  rsife  cos  0')2  +  r^ide  sin2  0'  +  (R  cos  9  -  z')2 


=  t-R/c- 


rside  sin  6  cos  0'  +  z'  cos  6 


in  the  far  field  region. 


When  evaluating  Equation  18  the  1/R  dependence  is  suppressed.  To  calculate 

side 

E  (R,  t),  H ^(rside,  z',  V)  and  E2(R',  f)  are  required  for  times  t  =  nAf.  Since  E  and  H 

side 

on  the  surface  are  only  calculated  at  times  0,  At,  2At,  then  E  (R,  t)  cannot  be 

I 

evaluated  for  any  given  t.  Let  t  ■  be  defined  as 


rside  sin  9  cos  0  +  z  ■  cos  9 


t-„  =  n  Af 

jnin  1 


where  the  R/  c  time  delay  has  been  ignored.  Consider  the  first  integral 


2k 

ztop 


UR0  = 


J Jsside,  z\  f)  dz'  dp 


Using  sum  approximations  for  the  0  and  z  integrals  leads  to 

2k  jtop 

IH(Rtk)  =  ^4>  ^2  ^  sider  j>  tjmn ) 

m=0  j=jgp 

where  the  time  is  defined  using 

k  =  int  (tjmn/  Af) 

*k  =  kAt 

Therefore,  the  field  at  a  given  point  on  the  side  on  the  cylindrical  surface 
z  =  Zj,  0  =  0m,  at  time  t  =  n  At  contributes  to  the  far  field  at  a  given  R  in  the  kth 
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retarded  time  "bin"  corresponding  to  time  tp.  The  integral  is  evaluated  by  summing 
the  contributions  over  all  times  from  all  points  on  the  surface  for  each  f \  bin.  All  of 
the  required  integrals  were  evaluated  in  this  fashion. 

At  the  top  of  the  cylindrical  surface 

n  x  H(R',  t')  =  -  H^R1,  t ')  [sin  6  cos0'R  +  cos  6  cos  0'  0  +  sin0'  (j)] 

and  hence 


RxRxnx  H(R’,  f)  =  H  (R1,  f)  [cos  d  cos  0'  6  +  sin0'  <))] 

Also 


n  x  E(R',  f)  =  Er(R',  t')  [  -  sin  6  sin  0'  R  -  cos  6  sin  0!  0  +  cos  0'  <[)] 
which  leads  to 

Rxnx  E(R',  f)  =  -  Er(R',  t')  [cos  0'  0  +  cos  6  sin  q Y  4>] 

The  contribution  to  the  far  zone  electric  field  from  the  top  of  the  cylindrical  surface  is 
therefore 


top  ix0  f A  A  d  A 

E  <R'f)  =  ^R[  JRxRx^[AxH(R \t')]dS' 

top 

'Vo  jRx|i[nxE(RVf ')]  rfS'] 

top 


4  xR 
In 

**  rside 

[  J  r '  H^(r'  >,  zf0p,  t')  dr'  [cos  0  cos  0'  9  +  sin  0'  0]  d(f>' 

0 
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2  K 

rrside 


1 


f  d  A  A 

r'  Er(R',  t')  dz'  [cos  0'  0  +  cos  0sin  <p'  <j>]  d(j)'  ] 


2  7t 

C  rside 


_w_ 

4  kR 


[cos  6 


C  d 

J  r'  ztop,  f)  dr'  cos  <p'  d<j>'  0 


2n 

rrside 


1 


r  d 

I  r'  ^Er(R',  t')  dr1  cos  0'  d<j>’  0  ]  (21) 


since  Er(R',  t')  is  an  even  function  of  0',  as  can  be  seen  from 

I  R  -  R'  | 
f  =  f- - L 

c 

V  (R  sin  6  -  r'  cos  01)2  +  r'2  sin2  <j>'  +  (R  cos  0  -  ztop)2 


=  f- 


=  t-R/c 


r'  sin  Q  cos  0'  +  Zf0p  cos  9 


As  for  Equation  20,  when  evaluating  Equation  21  the  1/R  dependence  is  suppressed. 


Let  f  ■  be  defined  as 


hmti 


I  I 

r;sin  0  cos  (j>m+  Zf0p  cos  6 
nAt_ - - - 


where  the  R/ c  time  delay  has  been  ignored.  Consider  the  first  integral 
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rside 

I(Rf)  =  J  r'  H  Jj'  i,  Zfgp,  t')  dr'  cos  0'  dtp') 

0 

Using  sum  approximations  for  the  <p  and  z  integrals  leads  to 

2k  hide 

I(Rtk)  =  A^Ar  X  X  COS  (mA<p)  jtopr  timn) 

m= 0  i  =  0 

where  the  time  tk  is  defined  using 

k  =  int  {tifitn/  Af) 

f*  =  fcA* 

Note  that  for  antennas  mounted  on  an  infinite  ground  plane,  the  fields  from  the 
antenna  image  below  the  ground  plane  must  be  included  in  the  surface  integration  to 
simulate  the  effect  of  the  ground  plane. 

From  the  expressions  derived  in  this  section,  the  electric  field  in  the  far  field  region  at 
an  angle  0  from  the  antenna  at  a  given  time  can  be  calculated  from  the  FDTD  field 
quantities  at  previous  times  on  a  cylindrical  surface  surrounding  the  antenna.  The  far 
field  calculation  is  central  to  obtaining  predictions  of  the  transmit  and  receive 
responses  of  the  antenna. 

2.6  Energy  Calculation 

To  calculate  the  efficiency  of  the  antenna  for  a  specified  excitation  voltage,  it  is 
necessary  to  integrate  the  power  radiated  from  the  antenna,  and  compare  it  to  the 
integrated  input  power.  By  checking  the  energy  reflected  from  the  antenna  back 
through  the  feed,  conservation  of  energy  can  be  used  to  test  the  validity  of  the 
calculation.  Following  [15],  the  net  power  transmitted  by  the  antenna,  E tmns>  is  given 
by 


oo 


where  V;„  is  the  applied  voltage  at  the  feed,  \refi  is  the  voltage  reflected  from  the 
antenna  back  through  the  feed,  and  Zcoax  is  the  impedance  of  the  coaxial  feed  line. 
The  efficiency  of  the  antenna  can  be  calculated  from 
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r  2  2 

Vin  -  Vrefl 

^coax 


dt 


n  = 


vL 


■dt 


-‘coax 


0 


The  power  radiated  through  the  same  cylindrical  surface  used  for  the  far  field 
calculation  was  determined  using 


rside 

pr«rf(0  =  J  2nr  E r(r,z,t)  H<p(r,z,t)  dr 
0 

over  the  top  of  the  cylindrical  surface  from  r  =  0  to  r  =  rsil}e,  and 

ztop 

side  r 

?rad  (0  =  -  J  2itr  E Z(r,z,t)  H0(r,z,f)  dz 

Z8P 

over  the  side  of  the  surface  from  the  ground  plane  at  z  =  zgp  to  z  =  zt0p.  The  electric 
fields  Er  and  Ez  were  linearly  interpolated  using  the  values  at  the  edge  of  the  grid  cells 
to  find  their  values  at  the  centre  of  the  cells  where  was  evaluated.  A  linear  time 
interpolation  was  also  needed  because  the  values  are  calculated  at  time  points 
(n+l/2)Af,  as  opposed  to  nAf  for  Er  and  Ez.  The  total  energy  radiated  by  the  antenna, 
E rad,  is  therefore 


r  top  side 

E rad  =  J  Prarf(0  +  Prarf  (0  dt 

oJ 

For  a  lossless  antenna,  E rad  should  equal  E tfans.  The  ratio  E^j/E^ns  was  calculated 
to  assess  the  validity  of  the  FDTD  calculation,  with  1  being  the  nominal  result 
provided  a  sufficiently  long  time  span  is  utilised. 

2.7  Coarse  Outer  Mesh 

In  order  for  the  absorbing  boundary  conditions  to  be  effective,  the  mesh  limits  were 
required  to  be  some  distance  from  the  antenna.  Since  the  fields  away  from  the 
antenna  are  well  behaved,  a  coarser  mesh  spacing  can  be  used  to  substantially  reduce 
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the  required  computational  overhead,  as  in  [12].  A  coarse  mesh  using  cells  of  size 
Nc Ar  x  NCAZ  were  used,  as  opposed  to  Ar  x  A2  for  the  fine  mesh  cells.  A  value  of 
Nc  =  3  was  generally  chosen.  The  interface  between  the  fine  and  coarse  meshes  is 
specified  by  the  user  through  iinner  and  jinner ■  The  alignment  of  the  meshes  is  shown 
in  Figure  2.7.  Note  that  the  upper  edges  of  the  lowest  cell  in  both  the  fine  and  coarse 
grids  was  aligned  with  z  =  0. 

To  calculate  E r(z,  jinner)  for  the  fine  mesh  along  the  top  of  the  cylindrical  interface, 
H<p(i/Nc,  jinner/ Nc  +  1)  from  the  coarse  mesh  and  H^z,  jinner  -  (Nc-l)/2)  from  the  fine 
mesh  were  first  used  for  values  of  i  which  were  multiples  of  Nc.  The  remainder  of  the 
fine  mesh  Er  values  were  obtained  from  these  values  using  linear  interpolation.  A 
similar  procedure  was  followed  for  Ez  along  the  side  of  the  cylindrical  interface  at 
i  =  Unner  -  (Nc  - 1)/2  - 1  for  the  fine  mesh.  The  same  time  steps  were  used  for  the  fine 
and  coarse  meshes.  Stability  for  the  coarse  mesh  follows  automatically  if  it  exists  for 
the  fine  mesh. 
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Figure  2.7:  Fine/coarse  mesh  interface  detail 
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2.8  Summary 

In  this  section,  the  techniques  used  for  the  rotationally  symmetric  FDTD  calculation 
were  discussed.  While  based  on  the  simple  Yee  algorithm,  the  contour  path  technique 
was  used  to  accommodate  arbitrarily  shaped  boundaries  without  resorting  to  the 
staircase  approximation.  While  this  approach  was  adopted  in  a  limited  sense  by 
Maloney  and  Smith  [21],  this  is  the  first  time  to  the  authors'  knowledge  that  this 
technique  has  been  applied  generally  in  cylindrical  geometry. 

A  simple  absorbing  boundary  condition  was  used  to  terminate  the  outer  limits  of  the 
mesh.  The  electric  field  in  the  far  field  region  was  calculated  by  integrating  the 
retarded  fields  over  a  cylindrical  surface  enclosing  the  antenna.  The  efficiency  of  the 
antenna  could  be  found  and  conservation  of  energy  verified  by  integrating  the 
radiated  power  over  the  same  cylindrical  surface  used  for  the  far  field  calculation.  To 
reduce  the  computational  effort  required  to  calculate  the  fields  sufficiently  far  from 
the  antenna  for  the  absorbing  boundary  conditions  to  be  valid,  a  coarse  mesh  was 
used  to  surround  the  inner  fine  mesh  covering  the  antenna. 

This  FDTD  calculation  was  implemented  in  the  program  fdtd_cyl.  In  the  next 
section,  results  obtained  from  this  analysis  for  spherically  capped  conical  monopoles 
are  compared  with  results  from  an  analytic  treatment.  Theoretical  predictions  for 
both  conical  and  cylindrical  monopoles  are  compared  to  experimental  data. 


3.  Conical  Monopole  D-dot  Sensor  Analysis 

3.1  Introduction 

A  small  conical  monopole  antenna  on  a  ground  plane  which  is  used  as  a  D-dot  sensor 
for  transient  electromagnetics  research  was  analysed,  with  parameters  a  =  0.625  mm, 
b  =  2.1  mm  for  the  coaxial  feed,  h  =  9  mm,  h'  =  0,  and  a  =  47°.  The  dielectric 
constant  in  the  coax  was  calculated  so  that  the  impedance  of  the  feed  line  was  50  Q, 
as  required.  This  meant  that  only  the  physical  dimensions  of  the  feed  line,  a,  and  b, 
needed  to  be  determined.  The  angle  a  corresponds  to  a  50  Q  characteristic  impedance 
for  the  monopole.  These  sensors  were  initially  described  by  Parkes  and  Smith  [22]. 

The  geometry  used  for  the  conical  monopole  is  shown  in  Figure  3.1.  An  optional 
hemispherical  cap  can  be  added  to  the  top  surface  of  the  cone.  The  cap  was 
modelled  using  the  simple  staircasing  technique.  This  was  found  to  be  sufficiently 
accurate  since  the  cap  was  sufficiently  far  from  the  feed  region.  The  contour  path 
technique  was  used  to  model  the  feed  and  the  conical  section  of  the  antenna.  Mesh 
cell  dimensions  were  Ar  =  Az  =  0.31  mm  for  the  fine  mesh,  with  iinner  =  200, 
jinner  ~  300,  and  imax  =  jmax  =  600.  The  number  of  fine  mesh  cells  per  coarse  mesh 
cell,  NCr  was  set  to  3. 
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Initially,  a  Gaussian  excitation  voltage  was  used  as  follows 

f-f2 1 

^inc(0  —  exP  2 

W 

with  Xp  related  to  the  FWHM  pulse  width  though  tw  =  2.355  xp. 

Q 


Figure  3.1:  Conical  monopole  geometry 

In  Appendix  A.l,  the  derivation  of  a  frequency  domain  mode  matching  treatment 
[25-27]  is  presented.  This  previously  developed  technique  was  used  to  provide 
independent  verification  of  the  FDTD  method.  In  addition,  the  calibration  of  the 
conical  D-dot  sensor  relies  upon  various  simplifying  assumptions  discussed  in  detail 
in  Appendix  A.l.  These  assumptions  are  used  to  derive  simple  results  from  the  full 
analytic  theory  which  relate  the  incident  electric  field  to  the  measured  output  voltage 
on  the  sensor  provided  the  sensor  is  small  compared  to  the  smallest  significant 
wavelength  of  the  incident  field. 

3.2  Reflected  Voltage 

From  Equation  A.17  in  Appendix  A.l,  the  input  impedance  of  an  ideal  electrically 
small  conical  monopole  can  be  approximated  by  a  capacitance  Ca,  where 
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Ca 


/z(4rcZc  +  3  £o  cos20o) 
4  nc2Z^ 


with  C,q  the  impedance  of  free  space,  and  Zc  the  characteristic  impedance  of  the 
monopole. 

The  frequency  domain  input  impedance  for  the  monpole,  Z^if),  was  calculated  from 
both  the  full  analytic  treatment  and  the  capacitance  of  an  ideal  monopole.  The  time 
domain  reflected  voltage,  Vrefi(t),  from  the  antenna  back  down  the  feed  line  was 
found  for  a  given  incident  voltage  Vinc(t)  using 

4*oo 

V„0)  -  J  Vinc(fl  expqatffl  if 


where 


4-oo 


/  ^ znc  (0  exp(-j2jc/f)  dt 


The  reflected  voltage  from  the  full  and  ideal  analytic  calculations  are  compared  with 
the  prediction  from  the  FDTD  method  in  Figures  3.2,  3.3,  and  3.4  for  values  of 
Xp  =  250,  100,  and  50  ps,  respectively.  For  all  three  values  of  xp,  the  agreement 
between  the  FDTD  and  full  analytic  results  is  excellent.  This  in  turn  implies  that  the 
mode  matching  calculation  correctly  models  the  input  impedance  of  the  small  conical 
monopole. 

The  agreement  between  the  ideal  results  and  the  other  two  methods  becomes 
progressively  worse  as  the  pulse  width  decreases,  as  expected,  since  the  antenna 
impedance  diverges  from  a  simple  capacitance.  Since  these  D-dot  sensors  are 
required  to  work  with  pulse  widths  as  low  as  100  ps  (xp  =  42.46  ps),  this  discrepancy 
in  the  ideal  reflected  voltage  cannot  be  disregarded  when  calculating  the  response  of 
the  sensor. 
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Time  (ps) 


Figure  3.2:  Reflected  voltage  from  ideal,  full  analytic,  and  FDTD  calculations  for  rp  =  50  ps 


Time  (ps) 


Figure  3.3:  Reflected  voltage  from  ideal,  full  analytic,  and  FDTD  calculations  for  Tp  =  100  ps 
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Figure  3.4:  Reflected  voltage  from  ideal,  full  analytic,  and  FDTD  calculations  for  xp  =  250  ps 


3.3  Transmit  Response 

In  the  frequency  domain,  the  transmission  transfer  function  of  an  antenna  is  given  by 


T(/) 


Eradjf) 


.CO f  he{f) 

]  cr  ZA(f)  +  Zg 


where  Vg(f)  is  the  open  circuit  source  voltage,  Eradif)  is  the  radiated  far  field,  ZA(f)  is 
the  antenna  impedance,  Zg  is  the  50  12  source  impedance,  and  he(J)  is  the  effective 
height  of  the  antenna.  The  open  circuit  voltage  Vg  is  twice  the  input  voltage  on  the 
coaxial  line  feeding  the  antenna  if  the  source,  feed  line,  and  antenna  characteristic 
impedance  are  matched.  The  parameters  ZA(f),  and  he(f)  could  be  calculated  for  the 
conical  monopole  using  Equations  A. 14,  A. 15,  and  A.18  from  Appendix  A.l.  Simple 
approximations  are  available  for  ZA(f),  and  he(f)  in  the  ideal  case  of  an  electrically 
small  monopole. 


As  discussed  in  A.6,  the  time  domain  radiated  electric  field  Erad(t)  on  the  ground 
plane  (9  =  90°)  in  the  far  field  region  can  be  calculated  for  the  full  analytic  and  ideal 
cases.  These  results  were  compared  to  the  predictions  from  the  FDTD  analysis  for 
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Xp  =  50,  100,  and  250  ps  in  Figures  3.5,  3.6,  and  3.7,  respectively.  The  excitation 
voltage  was  a  1  V  Gaussian  pulse.  The  1/R  dependence  has  been  suppressed,  where 
R  is  the  distance  to  the  far  field  point.  To  first  order,  the  simple  ideal  time  domain 
field  for  each  Xp  calculated  using  the  further  assumption  discussed  in  Section  A. 6 
that  the  monopole  impedance  Z&  »  Zcoax.  This  approximation  leads  to  the  simple 
expression 


3^q/i2cos  8q  d2Vs(t) 

rad\  )  47tZcrc2  dt2 

for  the  first  order  time  domain  radiated  field. 

From  Figures  3. 5-3. 7,  it  is  clear  that  good  agreement  exists  between  the  predictions 
for  the  far  field  derived  from  the  FDTD  and  analytic  techniques.  The  variation  of  the 
electric  field  Eq  in  the  far  field  region  predicted  by  both  methods  is 

Eq(0)  =  Eq(0)  sin  0 

as  expected  for  a  small  monopole.  The  ideal  analytic  calculation  using  small 
monopole  approximations  for  the  antenna  impedance  and  effective  height  works 
satisfactorily,  but  begins  to  break  down  at  shorter  pulse  widths  as  the  small 
monopole  approximations  break  down. 


Figure  3.5:  FDTD ,  full  and  ideal  analytic,  and  first  order  far  field  comparison  for  xp  =  50  ps 
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Time  (ps) 


Figure  3.6:  FDTD,  full  and  ideal  analytic,  and  first  order  far  field  comparison 
for  Tp  =  100  ps 


Time  (ps) 


Figure  3.7:  FDTD,  full  and  ideal  analytic,  and  first  order  far  field  comparison 
for  Tp  =  250  ps 
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3.4  Receive  Response 

The  receive  voltage  on  the  D-dot  sensor  can  be  calculated  from  the  transmit  response 
using  reciprocity,  since  the  frequency  domain  receive  transfer  function  S (/)  for  an 
antenna  can  be  expressed  as 


S  if) 


VilfL 

E  incif) 


Kif) 

zA(f) +  Zl 


where  Vi  is  the  receive  voltage  on  the  antenna,  ZA  is  the  antenna  impedance,  and  Zi 
is  the  50  Q  load  impedance.  For  Zi  =  Zg  from  section  3.3,  T (/)  and  S (f)  are  related  by 


s (/)  "  i~T (/) 


Therefore,  the  probe  voltage  Vi  for  an  incident  field  Einc  can  be  calculated  for  a 
receive  antenna  by  considering  the  same  antenna  in  a  transmit  mode.  The  radiated 
field  Erad  can  be  calculated  from  a  source  voltage  Vg  which  is  related  to  E,> lc  for  the 
receive  antenna  by 


V  (t)  =  — 

2  Tier  dt 


Under  these  circumstances,  the  receive  voltage  Vi  will  be  equal  to  the  calculated  Eraci 
for  the  transmit  antenna. 


For  receive  monopoles  which  are  small  compared  to  the  smallest  significant 
wavelength  in  the  incident  field,  the  voltage  developed  across  the  output  of  the 
antenna,  Vi  is  related  to  the  incident  field  E,„c  to  first  order  by 


VL 


3E2cos  9o  dEjnc(t) 
2c  dt 


if  the  angle  of  arrival  of  the  incident  field  is  90°,  ie.  parallel  to  the  ground  plane.  This 
equation  and  the  underlying  assumptions  required  to  derive  it  are  described  in  detail 
in  Appendix  A.l.  For  a  9  mm  conical  monopole,  the  constant  factor  is  0.2764  ps/ m, 
as  in  [22]. 

Because  of  the  reciprocity  relationship,  the  behaviour  of  the  receive  responses  follows 
that  of  the  transmit  responses,  as  shown  in  Figures  3.8,  3.9,  and  3.10  for  the  specified 
values  of  ip.  The  incident  field  is  a  vertically  polarised  Gaussian  pulse  with  a  peak 
amplitude  of  1  V/m. 
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As  for  the  transmit  response,  the  agreement  between  the  full  analytic  and  FDTD 
predictions  is  very  good.  The  ideal  analytic  calculation  using  small  monopole 
approximations  for  the  antenna  impedance  and  effective  height  becomes  somewhat 
inaccurate  for  shorter  pulses,  as  the  small  monopole  approximations  breakdown.  For 
T p  =  50  ps,  the  first  order  result  shows  the  worst  agreement  with  the  alternative 
techniques.  The  breakdown  of  the  first  order  response  relates  to  the  fact  that  the 
approximation  Za  »  Zcoax  is  no  longer  true  for  the  higher  frequency  components  of 
the  shorter  incident  pulse.  Since  the  approximation  Za  >>  Zcoax  is  fundamental  in 
determining  the  simple  D-dot  sensor  calibration  factor  used  in  [22],  this  can  lead  to  a 
significant  calibration  error  for  shorter  pulses.  Note  that  the  apparently  good 
agreement  between  the  FDTD  and  first  order  analytic  techniques  at  Xp  =  100  ps  is 
merely  fortuitous,  since  the  first  order  response  is  obtained  from  the  ideal  analytic 
response  with  further  approximation. 


Figure  3.8:  FDTD,  full  and  ideal  analytic,  and  first  order  receive  voltage  for  Tp  =  50  ps 
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Figure  3.9:  FDTD,  full  and  ideal  analytic,  and  first  order  receive  voltage  for  rp  = 
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Figure  3.10:  FDTD,  full  and  ideal  analytic,  and  first  order  receive  voltage  for  t 
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3.5  Comparison  of  FDTD  D-dot  Response  with  Experiment 

An  experiment  was  performed  to  verify  the  accuracy  of  the  FDTD  predictions  for  the 
9mm  conical  monopole  D-dot  sensor  receive  response.  A  D-dot  sensor  was  mounted 
on  a  ground  plane  31  cm  from  a  cylindrical  monopole  with  a  radius  of  6.25  mm  and 
a  length  of  40  cm.  The  field  at  the  D-dot  site  was  calculated  using  the  FDTD 
technique.  This  field  was  then  used  to  calculate  the  output  voltage  from  the  9  mm 
conical  D-dot  sensor.  The  predicted  D-dot  response  was  therefore  dependent  on  the 
accurate  FDTD  modelling  of  both  the  cylindrical  monopole  and  the  conical  sensor. 

A  Kentech  HMP  1  pulse  generator  was  used  to  drive  the  cylindrical  monopole  with  a 
3.5  kV,  2  ns  pulse  with  a  100  ps  risetime.  A  Barth  142-NMFP-26  attenuator  was 
used  to  reduce  the  Kentech  output  pulse  by  26  dB.  The  D-dot  sensor  output  voltage 
and  the  Kentech  output  pulse  were  measured  using  a  Tektronix  TDS  820  6  GHz 
sampling  oscilloscope.  Data  acquisition  was  performed  using  a  486  PC  with  a  GPIB 
interface  and  LabWindows  CVI  software. 

The  measured  and  predicted  conical  D-dot  sensor  output  voltages  are  shown  in 
Figure  3.11.  Note  that  some  high  frequency  components  of  the  applied  voltage  pulse 
have  been  attenuated  in  the  measured  data.  Note  that  the  FDTD  calculation  did  not 
include  lossy  dielectrics  or  conductors.  The  agreement  between  the  FDTD  results  and 
the  measured  data  is  satisfactory,  especially  since  two  antennas  have  been  modelled 
in  the  process  of  this  comparison. 


Time  (ps) 


Figure  3.11:  Measured  and  predicted  conical  D-dot  sensor  output  voltage 
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The  electric  field  from  the  cylindrical  monopole  as  calculated  by  the  FDTD  technique 
is  shown  with  the  reconstructed  electric  field  from  the  conical  D-dot  output  voltage. 
The  first  order  analytic  expression  in  Appendix  A.6  derived  from  the  mode  matching 
analysis  is  used 


t 

t 

=  3.620  x  10-12  J  VL(t')  dt' 

0 

In  Figure  3.12,  note  that  the  reconstructed  electric  field  underestimated  the  calculated 
field.  This  is  consistent  with  the  discrepancy  in  Figure  3.11.  In  addition.  Figure  3.8 
suggests  that  for  faster  pulses,  the  first  order  analytic  approximation  will  lead  to  an 
underestimate  of  the  measured  field.  Nevertheless,  the  agreement  is  considered 
sufficiently  good  for  practical  purposes. 


Figure  3.12:  Reconstructed  and  predicted  incident  field 
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4.  Conclusions 

A  finite  difference  time  domain  (FDTD)  technique  has  been  developed  to  analyse 
rotationally  symmetric  antennas  such  as  biconical  antennas  of  arbitrary  cross  section. 
The  standard  Yee  formalism  was  implemented,  with  a  cell  deformation  technique 
used  to  accurately  model  the  feed  region.  This  was  necessary  to  accurately  predict 
the  high  frequency  response  of  the  antenna.  A  coarse  grid  surrounding  an  inner  fine 
grid  was  used  to  reduce  computational  effort  for  outer  regions  away  from  the 
antenna.  Although  the  FDTD  technique  calculates  only  the  antenna  transmit 
response  directly,  the  receive  response  is  readily  calculated  using  reciprocity. 

A  9  mm  conical  monopole  D-dot  sensor  on  an  infinite  ground  plane  was  analysed 
using  both  the  FDTD  technique  and  an  analytic  technique  based  on  mode  matching. 
The  agreement  between  the  two  methods  was  excellent,  with  the  predicted  9  mm 
conical  D-dot  sensor  response  verified  by  experimental  results.  The  accuracy  of 
approximate  techniques  based  on  the  analytic  mode  matching  approach  was 
investigated,  with  some  discrepancies  with  the  FDTD  technique  noted  for  shorter 
pulses.  Now  that  its  accuracy  has  been  verified  by  the  FDTD  method,  the  full  mode 
matching  analysis  will  be  used  to  more  accurately  reconstruct  incident  electric  fields 
from  9  mm  conical  D-dot  measurements  for  shorter  pulses. 

Future  work  will  be  directed  toward  interfacing  the  routine  fdtd_cyl  to  CAD  software 
for  more  convenient  design  and  analysis  of  rotationally  symmetric  antennas  with 
arbitrary  cross  section. 
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Appendix  A 


Spherically  Capped  Conical  Monopole  Analysis 

The  frequency  domain  response  of  a  spherically  capped  conical  monopole  antenna  on 
an  infinite  ground  plane  as  shown  in  Figure  A.l  can  be  calculated  analytically.  The 
response  of  a  biconical  antenna  with  the  same  dimensions  can  be  determined  trivially 
from  the  monopole  analysis.  The  fields  within  the  spherical  boundary  of  the  antenna 
are  matched  to  the  fields  outside  the  antenna,  which  is  expressed  in  terms  of  spherical 
eigenmodes.  This  allows  the  reflection  coefficient  for  the  dominant  TEM  mode  at  the 
aperture  to  be  calculated,  which  in  turn  yields  the  antenna  driving  point  impedance 
and  the  far  field  radiation  pattern.  Because  of  the  requirement  for  a  spherical 
boundary  for  field  matching,  this  analysis  is  restricted  to  conical  monopole  antennas 
with  spherical  caps.  An  assumption  is  made  that  there  is  no  reflection  from  the 
transition  between  the  coax  and  the  antenna  feed  point. 


Figure  A.l:  Conical  monopole  geometry 
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A.l  Spherical  Eigenmodes 

The  spherical  eigenmodes  inside  and  outside  the  antenna  can  be  determined  from  the 
electric  Hertzian  potential,  since  only  TM  modes  are  excited.  The  electric  and 
magnetic  fields  are  written  in  terms  of  the  Hertzian  potential  ne  as 

H=  j<D£Vxne  and  E  =  ot^poETIe  +  V(V.ne) 

Since  Hr  =  0  for  TM  modes,  ne  has  only  a  radial  component,  so  that 

rie  =  nr  ur  =  (nj/)ur 

where  y  is  a  solution  of  the  scalar  Helmholtz  equation 

V2\jr  +  /c2v|f(r,0)  =  0  where  k 2  =  co2poe  [23]. 

Therefore 


1 A  (Jan  .  1  a 


^d~r{r2dr  J  +  if(sin  0  30  J +  k =  0 

This  equation  is  separable,  so  using  ijr(r,0)  =  R(r)  d(0)  results  in  two  equations 
„  d2R  dR 


dr2 


+  2r^+(fc2r2-p2)R(r)  =  0 


d2d  cos  0  dd  ~ 

dP  +  iiTede  +  P  ®*9)  *  0 

where  p  is  the  constant  of  separation.  Substituting  for  x  =  cos  0  in  (2)  yields 

d2d  dd 

(1'X)d^'2xd^+V(V+1)^X)  =  ° 


(A.l) 

(A.2) 


(A.3) 


where  v(v+l)  =  p2. 

The  solutions  to  Equation  A.3  are  the  Legendre  functions  of  the  first  and  second 
kinds,  Pv(cos  0)  and  Qv(cos  0).  These  functions  are  [24] 


Pv(*) 


=  Vrc 


Fj-y,  v+!;  1 ;  -t2)  _ 
1  v  v 

r(2-2)r(i+2) 


F(2"2/  ^  +  ?'  ?'  x 


2'  2’ 


1  v  v 

r(2+2)r(-9) 


42 


DSTO-RR-0019 


QvW  - 


Vtt 


1  v.  ,(V+l)jT  V  1  V  1  V  V7I  1  V  „  V  3  " 

F(2+2)  cos(  2  )  F(- 2 /  2  +  2'  2  '  X  ^  F(l+2)  cos(  2  )  F(2_2/ 1  +  2'  2'  *  ) 

- _ - +  x - -  — - 


2  r(i+|) 


r<H> 


where  the  hypergeometric  function  is 


F  (a,  b;  c ;  x) 


oo 

r(c)  V  T(a  +  n)r(b  +  n )  xn 
mm  ^  r (c  +  n)  n\ 


The  solutions  to  Equation  A.l  are  the  spherical  Bessel  functions  j \{kr)  and  y v(fcr), 
where 


jv(z)  and  yv(z)  can  be  combined  into  spherical  Hankel  functions  of  the  first  and 
second  kind  to  represent  waves  travelling  inward  and  outward,  respectively,  relative 
to  the  origin. 


=  jv(z)  +  jyv(z) 


inward  travelling  wave 


^2)(z)  =  jv(z)-iyv(z) 


outward  travelling  wave 


Since  the  field  outside  the  spherical  boundary  of  the  conical  monopole  must  be  finite 
at  0  =  0  and  rt,  v  can  only  take  integer  values,  n.  The  external  field  can  also  only 
consist  of  outgoing  spherical  waves,  the  scalar  function  \|/out  is  given  by 


Vout 


-  - <2) 


“  icoe 
«  =  o  J 


An  \  (kr)  Pn(cos  0) 
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So 


n 


r  =  X  J~rAn\if\kr)  P„(cos  9) 


„ =o 


where  A  are  arbitrary  and  the  l/j(oe  factor  has  been  included  for  the  sake  of 
consistency  with  the  derivation  in  [25]. 

a 

Since  =  0,  and  the  TM  mode  magnetic  field  has  only  a  (j)  component 

i  anr 

H  =  H<t,  =  -  jwe  -  -g§- 


n  =  0 


(2)  d 

Anh«  (^r)  d0  ^”(COS  ®) 


(2) 

An  h„  ( kr )  sin  0  Pn(cos  0) 


M  =  0 


The  electric  field  outside  the  conical  monopole  is  readily  determined  from  the 
magnetic  field  using 


Consider 


Since 


E  = 


: — V  xH 

jcoe 


E  r  = 


1 


jfcr  sin  0 


M-0  a 
e“  a©  (n<t> sin  e) 


d  f  dP„(cos0)" 

d0  vsin  0  d0  J 


=  sin  0 


d2P„(cos  0)  dPn(cos  0) 
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=  -sin  0  cos  0  P„(cos  0)  +  sin3  0  P„(cos  0)  -  sin  0  cos  0  P„(cos  0) 

f  11 

=  -  2  sin  0  cos  0  P„(cos  0)  +  sin3  0  P„(cos  0) 

Er  can  be  calculated  from 

rj  _  .^[E_L 

Er  “  \  e  kr  sir 


1  sin  0 


(2) 


«  =  0 


A„  h,,  (kr)  (2  sin  0  cos  0  P„(cos  0)  -  sin3  0  P„(cos  0)) 


(2) 


fao  V  Ki  M 

=  -  j  \  T  ^  A«  fo-'-  n(n+l)  sin  0  P„(cos  0) 


n  =  0 


smce 


(i-x2)  P„(x)  -  2xPm(x)  +  n(n+ 1)  P„(x)  =  0 
Similarly  Eq  can  be  determined  from 


Ee  = 


with 


j  kr  y  e  dr 


=  hf)(kr)  +  fcrhf)(*r) 

=  fcrh^^fcrJ-nhfVr) 


Therefore 


Ee  =  j  sin  0  P„(cos  0)  (kr  hf\(kr)  -  n  hf\kr)) 


n  —  0 


oo 

=  j  X  A„  sin  0  P„(cos  0)  (lV|(fcr)  -  ~  h{^(kr)) 


n  =  0 
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Since  the  field  inside  the  spherical  boundary  of  the  conical  monopole  is  not  required  to 
be  finite  at  9  =  0  and  n,  and  consists  of  both  ingoing  and  outgoing  spherical  waves, 
the  scalar  function  vyin  is  given  by 

oo 

Vin  =  X  (Av  hj,  \kr)  +  Bv  \\\kr))  (Py(cos  0)  +  avQv(cos  0)) 

Although  the  mode  spectrum  is  discrete,  note  that  v  is  not  restricted  to  integer  values. 
The  coefficients  Ay  refer  to  modes  propagating  outward  from  the  origin,  while  the  By 

coefficients  relate  to  modes  reflected  from  the  aperture  and  travelling  inward  towards 
the  origin.  Using  the  scalar  function  \|/in  and  proceeding  as  for  the  fields  outside  the 
spherical  boundary  of  the  conical  monopole  yields  the  fields  inside  the  interface  as 
follows 

OO 

V  (2)  (1)  d  d 

n<i>  =  -  Zj  ( W  (kr)  +  BV\  (kr))  (depv(cos e) +  av d^2v(cos  0)) 

v  =  0 

E  r  =  'j 

00  (2)  (1) 

Y  K  (kr)  K  \kr) 

Zu  (Ay  kr  +  Bv  . )  v(v+l)  sm  0  (Py(cos  0)  +  ayQy(cos  0)) 

v  =  0 


OO 

V  (2)  V  (2)  (1)  v  (1) 

L  (Ay(hy.j(fcr)  -  Vr  hi  (kr))  +  By(hi/1(fcr)  -  ^  ^  \kr))) 

v  =  0 

(|,Pv(cos  0)  +  av  ^2v(cos  0)) 

Values  of  v  corresponding  to  TM  mode  solutions  are  determined  from  the  boundary 

condition  on  the  side  of  the  conical  conductors,  Er  =  0  at  0  =  0Q  and  0  =  n  -  0Q. 
Therefore 

v(v+l)  sin  0  (Py(cos  0O)  +  ayQy(cos  0Q))  =0  (A. 4) 

and 
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v(v+l)  (Pv(-cos  e0)  +  avQv(-cos  0O))  =  0 


with  symmetry  requiring 


(A.5) 


v(v+l)  (Pv(0)  +  avQv(0))  =  0  (A. 6) 


These  conditions  are  trivially  satisfied  for  the  lowest  order  TM  mode,  the  TEM  mode 

dP0 

(v  =  0).  There  is  no  contribution  to  the  fields  from  the  P0(cos  0)  term  since  =  0. 
Therefore  the  TEM  fields  inside  the  conical  monopole  are 


TEM  (2)  (1)  d 

H«t,  =  -  (A0  hfl  (kr)  +  Bo  ho  (kr))  o(cos  e) 

=  A°  (exP(l*r)  '  P  exP0'fcr)) 

TEM  _  /llo"  (2)  (1)  d 

E0  =  - j V  e  (AohJiV)  +  B0hV(fcr))  ^Q0(cos  0) 

...  ;A/^a  exp(-jfer)  +  p  expQTcr) 

-J  \  e  A0  fcr sin 0 


(A  .7) 


(A.8) 


where  the  a0  factor  has  been  absorbed  into  A0  and  B()/  p  =  B0/  A0  is  the  reflection 
coefficient,  and  the  following  expressions  have  been  used 

_  ,  „v  1 ,  fl  +  cos  0" 

Qo(cos0)  =  jHi-cose 


-  ^  exp(j \kr) 
to  exp(-jicr) 

The  evaluation  of  higher  order  TM  modes  will  require  the  calculation  of  the  mode 
number  n.  This  is  performed  by  solving  the  boundary  condition  Equations  A.4,  A.5, 
and  A.6. 


-|qo(cos0)  =  ^ 

h o\kr)  =  -  j  h!J}(/cr)  = 

,(2)  x  .  A2) 
ho  M  =  J  h  j  (kr)  = 
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A.2  Aperture  Field  Matching 

In  practical  conical  monopole  antennas,  all  modes  inside  the  antenna  except  for  the 
TEM  mode  are  cutoff.  Since  this  mode  cannot  exist  outside  the  antenna,  power  is 
radiated  by  higher  order  TM  modes  which  are  excited  by  the  TEM  wave  at  the 
aperture.  The  scattering  coefficients  which  determine  the  amplitudes  of  the  TM 
eigenmodes  outside  the  antenna  can  be  determined  by  matching  the  transverse 
component  of  the  electric  field,  Eq  at  the  aperture  (r  =  a).  Inside  the  antenna  the 
approximation  made  by  Pappas  and  King  [25,  26]  is  used  which  ignores  the 
contribution  from  all  higher  order  evanescent  modes  inside  the  spherical  interface. 
This  approximation  will  be  a  fundamental  limit  to  the  accuracy  of  this  technique 
unless  higher  order  modes  are  included.  The  electric  field  inside  the  antenna  is 
therefore  given  by 


Ee 


-  ■  j  Ao 


ho  exp(-jfca)  +  p  exp(jfca) 
e  ka  sin  0 


over  the  aperture  (90  <  9  <  7T-@o) 


=  0  over  spherical  cap 

where  p  is  the  reflection  coefficient  which  accounts  for  TEM  mode  reflection  at  the 
aperture  discontinuity.  The  field  outside  the  antenna  is 

oo 

Ee  =  -  j  X  A„  (h^(fcfl)  -  b^(ka))  ~^Pm(cos  0) 

n  =  0 


Equating  the  fields  inside  and  outside  the  aperture  leads  to 

OO 

A0  X"1  (2)  m  (2)  d 

la  sin  Q  (exp(-jfca)  +  P  exp(jfca))  =  ^  Am(h -  ^hw  (*«))^Pm(cos  0) 

m  =  0 

Following  the  method  of  Papas  and  King  [25,  26],  multiplying  both  sides  by  sin  0 
P„(cos  0)  and  integrating  over  the  spherical  interface  from  0  =  0  to  jt  yields 


7t-0o 

A0  (exp(-jfcfl)  +  p  expQka))  f  (^P„(cos  0))  d0 

eo 


2  A0  P„(cos  0O)  (exp(-j ka) 


+  p  exp  (pi)) 
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oo  K 

X  ka  A M(h21(*B)  -  Ya  ^\ka))  J  sin  e  ^P„(cos  0)“  Pm(cos  0)  d0 

m  =  0  0 


»  2teA„f^(h£»(H-£hiV)) 

where  the  following  orthogonality  relation  has  been  used  [27] 


n 


J 


d  d  2n(n+l) 

sin0d0P»(cos  0)d0P>"(cos  0>  d9  =  2n+l  5mn 


Because  of  the  antenna  symmetry,  only  odd  n  need  be  considered.  Therefore  the 
spherical  eigenmode  amplitudes  in  the  region  outside  the  interface  are  given  by 


=  A0^r(-w;p^r(Wti^;)P„(eo,80) 


,  (2)  ,,  ,  n  (2) 


(A.9) 


Using  the  substitution  for  A„  in  Equation  A.9,  the  electric  field  outside  the  antenna 
can  be  related  to  the  reflection  coefficient  p  and  the  amplitude  of  the  outward 
travelling  TEM  mode  in  the  antenna  A0  by 


T7  —  •  a  ^/Boexp(-j/cfl)  +  p  exp(\ka) 

Ee  -  ■  J  A0  \  e  ka 

A2) ,  n  (2) 

V  2n+l  Vl  (kr)-frKikr)  d 

Vl2) - P"(COS  0o)^P«(cos  0) 


n(n+ 1)  (2)  n  ,  (I 

nodd 


(A. 10) 


'(ka) 


where  the  summation  is  over  odd  n  only.  In  the  far  field  region  where  kr  — >  °°,  the 
term 

From  Equation  A.10,  the  electric  field  in  the  far  field  region  is  therefore 
r,  -a  -*  exp(-jfcr)  exp(-\ka)  +  p  exp(j  ka) 

Ee  --jA0^e  kr  ka 
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CO 


I 

n  odd 


2n+l 


n(n+ 1)  (2)  n  (2) 


Pn(cos  eo)^Pn(c°S  0) 


(A.ll) 


The  electric  field  can  be  expressed  in  terms  of  the  current  Iq  through  the  centre 
conductor  at  the  feed  point.  Consider  a  circular  patch  of  radius  rc  on  the  ground 
plane  which  encloses  the  centre  conductor  of  the  monopole.  Ampere's  law  states  that 


iH- 


dl  =  I0  +  Yt  II  eE AS 


1Lee-‘ 


For  the  circular  patch 


2ti  rc  2k 

J  rcHtyd<\>  =  I0  +  ^  J  J  e  r  Ee  dty  dr 


0  0  0 
Because  the  monopole  is  azimuthally  symmetric 


27c  0<t,  r  =  lo  +  F (d,dt)  2ke  r  Ee  dr 

oJ 


Substituting  for  and  Eq  from  Equation  A.  7  and  A.8  leads  to 
\2n  j27i 

k  A0  (exp(-jfcrc)  -  p  exp(j krc))  =  I0  +  A0  [exp(-jfcrc)  -  p  exp(j krc)  -  (1-p)] 

At  the  feed  point  where  rc  — »  0,  the  second  term  on  the  right  tends  to  0,  so  that  the 
amplitude  of  the  outward  TEM  electric  field  in  the  monopole  is  related  to  the  applied 
current  by 


A  -  -i- 
0  j2n(l-p) 


By  substituting  for  A0  in  Equation  A.ll,  the  electric  field  in  the  far  field  region  can  be 
expressed  in  terms  of  the  feed  current  Io  and  the  TEM  mode  reflection  coefficient  p  as 

^  /  ftp  Ip  exp(-jfcr)  exp(-j ka)  +  p  exp(\ka) 
e  V  £  2jr(l-p)  r  ka 
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2n+l 


n(n+l) ,  (2)  n  ,  (2 


Jn  (2)  p«(cos  e0)^Pw(cos  0)  (A.12) 


\ka) 


This  is  identical  to  Equation  8  in  [27].  According  to  Harrison  and  Williams  [27],  the 
effective  height  of  the  conical  monopole  is  defined  as 


E  e(/) 


khe  ®  (jCo/  2n)  (exp(-j kr)/r)  1 0(/) 

where  =  V po/e  is  the  impedance  of  free  space. 
Using 


(A.13) 


h?V)  -  -r'yw 

as  kr  — >  oo,  the  magnetic  field  in  the  far  field  region  is  therefore 

_  Ip  exp(-jicr)  exp(-jfca)  +  p  exp(jicg) 

^  “  "2tc(1-p)  r  ka 


2n+l 


n(n+l) ,  (2)  „  x  n  (: 
"0<w  ;hi4(fcfl)-^h; 


J„  (2)  Pn(cos  eo)dd0Pn(cos  0) 


» 


E0 

=  ^ 

Provided  the  TEM  mode  reflection  coefficient  p  can  be  determined,  the  fields  outside 
the  antenna  can  be  calculated.  To  calculate  p,  the  tangential  component  of  the 
magnetic  field,  H^,  is  matched  over  the  aperture 


j  V  (2)  d 

ka  sin  0  (exP('jfcfl)  -  P  exP(j  ka))  =  Zj  KK  (ka)  d0P«(cos  e) 

n  =  0 

d 

Multiplying  both  sides  by  sin  0  Q0(cos  0)  and  integrating  over  the  spherical 
interface  from  0  =  0O  to  n-0o  leads  to 
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7l-0o 

■*"  ( exp(-jfc« )  -  p  exp(jfcg))  J  d0  = 


ka 


-  2  fa  ( exp(-jka )  -  p  exp(j ka))  In  cot 


V2  J 


Or 


7T-0 


exp(-jfcg)  +  p  exp(j ka)  V 

bn  ^ 


o 


n  =  0 


2n+l  .  |  ,d 


ka  Zj  ^^P«(cose0)U^)  J  (dep»(cos  e)) de 

e0 


2n+l 
n(n+ 1) 


[P„(cos  0O)]2 


where 


Uka) 


J2) 

K  (ka) 

hi-l  (ka)-^^\ka) 


Therefore 


(exp(-jfca)  -  p  exp(j/ca)) 


oo 


=  -  (exp(-jfca)  +  p  exp(j ka)) 


Ko_  V 

2*Zc  ~ 


2n+l 
n(n+ 1) 


[p„(cos  e0)]2 


^(ka) 


where  Zc  is  the  characteristic  TEM  mode  impedance  of  an  infinitely  long  conical 
monopole  given  by 


Consider 


Co, 

Z,  =  —  In  cot! 
IK 


fc,  \ 

yo 

UJ 


exp(-j2fca)  -  p  =  -  y  (exp(-j21ca)  +  p) 


so  that 
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P 


exp(-j  2 ka) 


1+Y 

1-y 


=  exp(-j  2 ka) 


N 


1  +  i2^^  S |^P>»(c°»(eo))]2W^ 


ti  odd 

~w 


(A. 14) 


where  the  summation  is  over  n  =  1,  3,  5,....,  N. 

Now  that  the  TEM  mode  reflection  coefficient  has  been  determined,  the  fields  outside 
the  conical  monopole  can  be  calculated  using  the  modal  coefficients  A„.  This  result 
has  been  presented  in  [25-27],  The  impedance  of  the  finite  bicone  antenna  is  therefore 

z A(f)  =  zc  (A. 15) 


A.3  Power  Conservation 

Conservation  of  power  across  the  spherical  interface  separating  the  regions  inside  and 
outside  the  spherical  boundary  was  used  to  assess  the  influence  of  neglecting  higher 
order  modes  within  the  antenna.  Inside  the  antenna,  the  TEM  electric  and  magnetic 
fields  can  be  used  to  calculate  the  power  crossing  the  interface,  Pr,  at  r-  a  using 


71/2 

Prn  =  7t a2  |  Re{Ee(fl,e)  H0(fl,0)*}  sin  0  d0 
00 


2  «/  2 

l0  Co  f  1 

4tt(1  -  p)2  0--P  +  RefP  exP(j2fc«)  -  P*  exp(-j21cfl)))  J 

00 

Io  Co  _  /0ox 

4jc(1  _p)2  (1  -  P2  +  Re{p  expO'2 ka)  -  p*  exp(-j2fca)})  In  cot  (yj 


IqZc1-P2  +  Re{p  exp(j2fc«)  -  p*  exp(-j2/ca)} 
2  (1-P)2 
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Note  that  if  there  is  no  reflection  from  the  interface,  then  all  of  the  applied  power  is 
radiated.  The  applied  power  is  given  by 


1  r 


2 

Iq  zc 
2 


Using  Equation  A. 12,  the  power  crossing  a  hemispherical  boundary  surrounding  the 
antenna  at  r  — »  «>  is  given  by 


jc/2 


out 


n r2  J  Re{Ee(r,0)  sin  0  dQ 


Io  Co 


4tc(1  _p)2  C1  +  P2  +  Re(p  exp(j2 ka)  +  p*  exp(-]2ka)}) 

N  N 

V  y  (2m+l)  (2n+l)  _ -  f Pm(cos(0o))  Pw(cos(0o)) 

mmmd  \ 


* » -  m(”+iw’,+1)  [ta  .  „  hiV)]  [*.  Cm  - m  hS'ctaj] 

rc/2 


J 


sin  0  ^Pw(cos  0)“  PM(cos  0)  dQ 


Since 


n 


f  d  d  2n(n+l) 

J  sm  0  d0P"(cos  0)d0  Pm(cos  0)  dQ  =  -2yt-+l  •  8mn 

the  radiated  power  in  the  far  field  region  is  given  by 
out  Io  Co 

Pr  =  4ti(1  -p)2  C1  +  P  +  Re(P  exp(j2^«)  +  p*  exp(-j2/ca)}) 


N 


y  (2n+l) 

mmmU  1-7  /  vr  4- 1  \ 


(-1  )n+1  [Pn(cos(Q0))32 


rt  odd 


n(n+ 1) 


[ka  h|f.](/c«)  -  n  h^  \/ca)  jj 
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Power  conservation  across  the  spherical  interface  separating  the  inner  and  outer 
regions  of  the  monopole  was  verified  to  a  high  degree  of  accuracy  in  the  frequency 
domain  up  to  at  least  ka  =  100  using  N  =  59  for  the  summation  terms. 

A.4  Small  Conical  Monopole  Approximation 

For  electrically  short  monopoles  with  ka  «  1,  a  number  of  simplifying 
approximations  can  be  used  to  calculate  the  input  impedance.  Consider  the 
summation  term  from  Equation  A.14 


For  small  ka 


Hence 


N 


s  = 


.  Co 

]2nZ, 


'c  n  odd 


fp«(cos(0o))]2  C«(*»)) 


^n\ka) 

-  fa  h?W) 

ka 


S  - 


[p«(cos(eo))]2 


For  0o  =  47°  corresponding  to  a  characteristic  impedance  of  50  £2,  then  the  dominant 
term  will  correspond  to  n  =  1.  The  next  highest  contribution  will  be  from  the  n  =  3 
term  which  is  only  1.5  %  of  the  n  =  1  term.  Therefore  S  can  be  approximated  by 


S 


■l^cos2(0o) 


The  reflection  coefficient  p  is  given  by 


p  =  exp(-j  2 ka) 


1  "  i|^fcflCOs2(e°) 

3  Co 

L1 + >M:ka  cos2(0»)J 


55 


DSTO-RR-0019 


1  -  ]ka 


1  +  )ka 


3^o 

1  ~i4xzcka  cos  (0o) 

3  Co 

L1  +  ii^fcflcos2(H 


1 '  (3i=°  cos2(0o)  +  4jtZc) 

ka 

1  +  (3^°  cos2(0°)  +  4jtZc) 


4jiZc  1 

jfcfl  3^o  cos2(0o)  +  4ttZc  ^ 
4tcZc  1 

\ka  3^o  cos2(0q)  +  47rZc  +  3 


Za 


- 1 


Za 

Zc 


+  1 


and  so  the  antenna  impedance  Za  corresponds  to  a  capacitance  Ca 

1 


ZA(/C«)  =  7 


j2Jt/CA 


(A. 16) 


where 


CA 


fl(4rcZc  +  3  cos2(0o)) 

2 

47rcZc 


From  Equation  A.12  and  A.13,  the  effective  length  he(f)  is  given  by 


(A. 17) 


khe(f)  =  - 


exp(-jfca)  +  p  ex 


ka  (1-p) 


?Qka) 


2n+l 


i  n-1 


"  odd  ^  ^  \(ka )  -  Ya  h^2)(/cfl) 


(cos  0o)^P„(cos  0) 


(A. 18) 


where  the  summation  is  over  n  =  1,  3,  5 For  electrically  short  monopoles 
(ka  «  1),  then  only  the  n  =  1  term  need  be  considered  in  the  summation,  as  for  the 
antenna  impedance.  Using  the  approximation 
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(2)  1  (2)  j 

ho  (ka)-  Shj  (to)  -  l-(t0)3 


(to)= 


the  effective  height  for  short  conical  monopoles  can  be  calculated  from 

,  3  exp(-j ka)  +  p  exp(ifca)  - 

khe(f)  ~  )  2  k a  (l-p)  (fcfl)3  cos  0O  sin  9 

Using  the  approximation 

exp(-jfca)  +  p  exp(jfcg)  1  -  j  ka  +  p  (1  +  j  ka) 

l-p  “  l-p 


i±£  , 


l-p 


-  jka 


together  with  the  approximation  for  p  in  (16)  leads  to 

exp(-jfcfl)  +  p  exp(j/ca)  _ 

l-p  ka 


j4tcZ^  (3^°  cos2  00  +  47tZc) 
4tiZc 

j ka  (3^o  cos2  0o  +  4jiZc) 


-]ka 


which  yields  the  effective  height 


6nkaZccos  0o  sin  0 
khe W  =  4tcZc  +  3Co  cos20o 


(A. 19) 


A.5  Transfer  Functions 

Following  the  standard  definition  in  [27]  the  transmission  transfer  function  is  defined 
as 


m 


Mi 

v0(/) 


where  Ee  is  the  electric  field  in  the  far  field  of  the  monopole  and  Vo  is  the  driving 
voltage.  Vq  is  related  to  Iq  in  (12)  by 
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T  -  V° 

Io  ZA  +  ZL 


where  Za  is  the  antenna  impedance  and  Zl  is  the  impedance  of  the  feed  line.  Using 
Equation  A.18,  T(/)  is  given  by 


W 


■  Co  /  he(f) 

]  cr  ZA(j)  +  ZL 


exp(-j2  nfr/c) 


(A. 20) 


=  T 0(J)  exp(-j2 nfr/c) 


To (/)  is  clearly  related  to  T (/)  with  the  phase  delay  due  to  the  distance  between  the 
antenna  and  the  far  field  point  removed. 

The  reception  transfer  function  S (/)  is  defined  in  a  similar  fashion  as 


S  (/) 


Vl  if) 

Einc(/) 


he(f)  ZL(/) 
Za  (J)  +  zl 


(A. 21) 


where  Vl  is  the  load  voltage  across  the  antenna  and  Ejnc  is  the  incident  field.  The 
effective  height  can  be  calculated  from  Equation  A.18  and  the  antenna  impedance 
from  Equation  A.15  as  a  function  of  frequency. 


For  small  monopoles  with  ka  «  1  which  are  matched  so  that  Zc  =  Zl,  the 
approximate  expressions  for  effective  height  Equation  A.19  and  antenna  impedance 
Equation  A.17  can  be  used  in  Equation  A.20  and  Equation  A.21  to  yield 


ToW  -  (j2 ltf)l 


3^o R2  cos  9o  sin  0 
47cZcrc2 


(A. 22) 


and 


3a2  cos  0nsin  0 

S  if)  =  j2  nf - Yc -  (A-23) 

if  the  assumption  ZA  »  Zl  is  made.  This  will  be  the  most  limiting  assumption  for 
the  case  of  a  50  Q.  monopole  connected  to  a  50  H  line. 

The  real  and  imaginary  components  of  the  transmission  transfer  function  To  for  a 
50  Q  monopole  with  0o  =  46.98°  are  plotted  against  ka  in  Figure  A.2.  The  1/  r  term  in 
Equation  A.20  has  been  suppressed.  This  plot  is  in  agreement  with  Figure  13  from 
[27],  Figure  A.3  shows  the  real  and  imaginary  components  of  the  reception  transfer 
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function  S (/)  for  50  fl  conical  monopoles.  The  agreement  with  Figure  16  of  [27]  is 
very  good. 


ka 


Figure  A.2:  Transmission  transfer  function  of  matched  conical  monopole 


ka 


Figure  A3:  Reception  transfer  function  per  unit  antenna  length  for  matched  conical 
monopole 
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A.6  Time  Domain  Response 

For  transmit  antennas,  the  time  domain  far  field  can  be  calculated  from  the  applied 
voltage  and  the  transmission  transfer  function  using 


-f-oo 

e0(£  +  r/c)  =  |  V0 (/)  T 0(/)  exp(j2nft)  df 

-oo 

The  frequency  domain  applied  voltage  is  calculated  from  a  time  domain  applied 
voltage  using  the  Fourier  transform 

+  00 

Vo (/)  =  |  vo(0  exp(-j2nft)  dt 


For  a  Gaussian  excitation  pulse  with  a  maximum  of  1  V,  the  applied  voltage  vo(£)  is 
given  by 


vo(£)  =  exp(-f2/2x2) 

where  x  =  0.4246  tW/  with  tw  the  FWHM  for  the  pulse.  The  applied  voltage  in  the 
frequency  domain  is 


Vo (/)  =  V  2tt  x  exp(-27i:2x2/2) 

For  small  monopoles.  Equation  A.22  can  be  used  to  relate  the  excitation  voltage  to 
the  transmitted  field  in  the  time  domain  as 


ee(f  +  r/c) 


3^ofl2  cos  8o  d2v0(f) 
4nZcrc2  dt2 


The  radiated  electric  field  for  a  number  of  antennas  with  different  a/  tw  values  is 
shown  in  Figure  A.4  for  a  Gaussian  excitation  pulse.  As  for  Figure  A.2,  the  1/r  term 
in  the  transmission  transfer  function  has  been  suppressed.  For  large  antennas  as  in 
Figure  A.4  with  a/  tw  =  6.62,  the  radiated  pulse  is  a  faithful  reproduction  of  the 
excitation  pulse.  For  short  antennas  as  for  a/tw  =  0.0883,  the  radiated  pulse  is  close 
to  a  double  derivative  of  the  excitation  pulse.  This  is  the  result  expected  from  the 
small  monopole  analysis.  The  results  presented  in  Figure  A.4  closely  match  those  of 
Figure  5  in  [27], 
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—  a/tw=  0.0883  m/ns 

—  a/tw=  0.442  m/ns 

—  a/tw=  2.21  m/ns 

—  a/tw=  6.62  m/ns 


Figure  A  A:  Radiated  field  from  50  Q  conical  monopole  for  various  a/tw 


The  time  domain  load  voltage  in  reception  can  be  calculated  from  the  incident  field 
using  the  reception  transfer  function  as  follows 

+  oo 

vl(0  =  J  E inc(/)  S if)  exp(j27t/0  df 

-oo 

The  incident  electric  field  is  calculated  in  a  completely  analogous  way  to  the 
excitation  voltage  for  transmission. 

For  small  monopoles.  Equation  A.23  can  be  used  to  relate  the  incident  field  to  the 
load  voltage  developed  across  the  antenna  in  the  time  domain  as 

3a2cos  6o  iteinc(t) 
vL(t) - £  rfT“ 

The  load  voltage  across  conical  monopoles  for  a  negative  Gaussian  incident  field  with 
a  peak  of  -1  V/m  are  shown  in  Figure  A.5.  Figure  A.5  shows  the  response  of  a  large 
monopole  with  a/  tw-  6.62  to  be  approximately  the  integral  of  the  incident  field.  For 
a  small  monopole,  the  response  is  proportional  to  the  derivative  of  the  incident  field, 
as  for  a/  tw  =  0.0883.  This  is  expected  from  the  small  monopole  analysis.  The  results 
shown  in  Figure  A.5  are  in  close  agreement  with  those  presented  in  Figure  10  of  [27]. 
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-  a^w~  0-0883  m/ns 

-  a/tw~  0.442  m/ns 

-  a/tw=  2.21  m/ns 

- a/f^  6.62  m/ns 


Figure  A. 5:  Voltage  across  50  Q  conical  monopole  receive  antenna  for  various  a/tw. 


A.7  Summary 

The  results  presented  in  this  appendix  show  conclusively  that  the  analytic  theory  of 
[25-27]  has  been  successfully  reproduced,  and  that  the  approximations  are  identified. 
The  most  serious  approximation  is  the  inclusion  of  only  the  TEM  mode  within  the 
antenna,  with  higher  order  evanescent  modes  neglected  completely  inside  the  antenna. 
Of  secondary  importance,  it  is  assumed  that  there  are  no  reflections  from  the  feed. 
This  approximation  is  significant  when  calculating  the  high  frequency  response  of  the 
antenna,  since  the  feed  geometry  is  the  limiting  factor.  This  is  particularly  important 
for  the  analysis  of  high  power  antennas  where  the  feed  is  generally  made  as  large  as 
possible  without  degrading  the  high  frequency  response  of  the  antenna  beyond  an 
acceptable  limit. 

The  simple  expressions  for  the  antenna  impedance  and  effective  height  of  small 
monopoles  with  ka  «  1  presented  in  [27]  were  derived.  These  simple 
approximations  have  been  used  to  calibrate  small  antennas  such  as  D-dot  sensors.  If 
the  further  approximation  that  the  antenna  impedance  is  much  larger  than  the  load 
impedance,  then  simple  time  domain  expressions  relating  the  incident  field  and  load 
voltage  for  reception  and  the  excitation  voltage  and  radiated  field  for  transmission 
could  be  derived.  This  last  approximation  is  the  most  limiting  in  terms  of  the 
minimum  pulse  widths  that  can  be  accurately  modelled  using  these  approximate 
time  domain  formulae. 
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